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ABSTRACT 


Much work has been done in recent years studying 
inelastic electron tunneling in metal-insulator-metal 
eunpel junctrons. sInelastic Electron Tunneling Spectros-— 
copy (IETS) has become a method for investigating various 
molecules within a tunnel junction's insulating layer. 

In particular, AX-insulator-Pb tunnel junctions, where 
the insulator is aluminum oxide deliberately doped with 
organic molecules, have become fairly popular. The IET 
spectra produced, however, are complex and decipherable 
primarily via knowledge of infrared and Raman spectroscopy. 
This thesis explores the undoped A&-AX oxide-Pb tunnel 
junction itself. By studying how the junction ages - 

or anneals - insight into the nature of the aluminum 
oxide insulator, and its interaction with the AX and Pb 
metal electrodes is gained. Also, by applying a voltage 
of either polarity (i.e. AX positive or negative) across 
the two electrodes during annealing, or by uSing a cover 
electrode metal other than Pb, one acquires additional 


knowledge about the junction. 
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CHAPTER 1 


INTRODUCTION 


1.1. Preamble 

This thesis sets out to describe some of the 
effects which occur in aluminum-aluminum oxide-lead 
(A£-oxide-Pb) tunnel junctions when subjected to thermal 
and voltage annealing. A-oxide-Pb tunnel junctions 
were chosen for several reasons: 1) they are the most 
commonly used type of junction for IETS - the oxide layer 
being doped with organic molecules, 2) ease of fabrica- 
tion, and 3) they yield "nice" spectra. 

The theory of tunneling, both elastic and inelastic, 
will be discussed. While the theory allows for the 
formulation of models which aid in the characterization 
of some of the changes occurring during annealing, no 
theory - from the tunneling point of view - has been 
advanced to really describe the time dependent effects 
which .do occur .in.tunnel. junctions, 

Junction preparation and the instrumentation used 
for measurement of junction properties will be described. 
Following this, data analysis will be dealt with. 

Finally, results of the annealing experiments will 
be presented along with models and conclusions derived 
from the experimental information. 


Suggestions for further work and appendices will 


conclude the thesis. 
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CHAPTER 2 


THEORETICAL CONSIDERATIONS 


Ze INneLOGUCL ion 

The concept of particle tunneling has been in 
existence almost as long as quantum mechanics itself. 
A particle, as represented by a wave function, may enter 
a classically forbidden region and if the potential 
barrier is sufficiently thin, has a certain probability 
om tunneling through. “That is, according to quantum 
mechanics, the particle may pass through the barrier 
without having enough energy to go over the top. 

Physically, this concept may be realized in the 
form of a tunnel junction in which electrons tunnel 
through an insulating barrier sandwiched between two 
metal electrodes. Fig. 2.1 shows the basic set-up of a 
tunneling experiment which allows one to measure the 
current-voltage characteristic of the junction. What 
happens if we attempt to pass electrons from one metal 
electrode to the other? Classically, one simply charges 
the electrodes as in a capacitor with no current flowing 
across the insulator. Quantum mechanically, however, one 
can get a current flow across this region provided that 
the insulating barrier is thin enough. This is electron 
tunneling. 

In order to understand the tunneling process we 


start by examining the electron energy level diagram for 
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Hal emer cel: 
This diagram shows the basic arrangement of a 
tunneling experiment. M is a metal electrode, 


I the insulator, V a voltmeter and A an ammeter. 


a normal metal-insulator-normal metal junction as shown 
Siig (aes Ch 9 ea ar EPl and En are the Fermi energies of 
metal 1 and metal 2 respectively while eV is the applied 
bias energy. Tunneling occurs from left to right (as 
denoted by the arrow) across a barrier of thickness s 
and barrier heights rT and 5 at the two electrodes. 

At T = OK all the electron states below Er. and 
E are filled, while those above are empty. The appli- 


F2 

cation of a voltage V across the junction (with the left 
hand metal biased negatively) causes the Fermi levels to 
become separated by an energy eV. Filled states of metal 
1 are now aligned with available empty states of metal 2 
and tunneling of electrons occurs from left to right. 
Most of this charge transfer occurs elastically (i.e. 
tunneling without electron energy change) although 
inelastic effects, which will be discussed later, do 
occur. As the voltage, V, is increased, more alignment 
of filled states with available states occurs and hence, 
the tunneling current is also increased. Thus, depending 


on the bias applied to the tunnel junction, one gets a 


certain current flow across the junction. 


2.2.) Blasticylunneling 

Most of the current flow across the tunnel junction 
is due to elastic tunneling in which the electron energy 
does not change when it passes through the barrier. The 


work of Harrison (59) which deals with tunneling from an 
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1D aRG a A 
The electron energy level diagram for a metal- 
insulator-metal junction is shown with and without 
an applied bias V. Here, s is the barrier 
thickness, 4 and d. are the barrier heights at 
the base and top electrodes respectively, and 
Erl and Ep? are the Fermi levels of the two metals. 


When a bias is applied, tunneling occurs in the 


direction indicated by the arrow. 


Metal 1 Insulator Metal 2 
(base) (top) 
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independent-particle point of view gives a theoretical 
calculation of the elastic tunneling current which will 
be briefly shown. Using Bardeen's work (10), the proba- 
bility per unit time for the transition of an electron 
from a state a in metal 1 to a state b in metal 2 on the 


other side of the potential barrier is given by 


ep RA 2 
epee (3 | poe peoet tots) G21) 


where MS is the matrix element for the transition, Pb is 
the density of states at b, and tf. and fy are the 
probabilities of occupation of states a and b 
respectively. 

Unless the transverse wave number ky is the same 
for the initial and final states the transition matrix 
element Mob vanishes and thus Py is a density of states 
for fixed k,. By summing over all states a of fixed ky 
to obtain Par summing over ky multiplying by 2 for spin 
and by e for electron charge, we obtain the total current 
fvom Jett to right. Subtrageing the jourrent an the 


opposite direction yields the current density 

.. pa A The 2 _ 

j= 4me 5 | M1? 030, (f4-f,) a2 (2.2) 
with the integral over energy being taken at fixed 


transverse wave number. 


In the evaluation of matrix elements, states are 
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constructed which are sinusoidal in a positive-energy 
region and drop exponentially in the adjacent negative- 
energy region. Assuming that the band structure is 
uniform, except near the transition region, and also that 
the band structure varies only slowly in the transition 
region allows us to make a WKB approximation. This 


results in the following expression for the current density 


oo s 
2e 


areal if tps | lk, |ax} (£,-£,) ak (2.3) 
t 0 0 

where s is the barrier width. We see that the density of 
states is absent from this expression; this is a conse- 
quence of the independent-particle model. 

By rewriting the is wavevector within the barrier 

, 1 +5 . 

region as ky = x [2m((x,V)-E,) ] and substituting the 
appropriate Fermi-Dirac distribution functions for ft. and 


it the current density may be written as 


et! 
“ s 
& 
ea <5 | exp{ - 2 | [2m($(x,V)-E,]* dx} x 
cere ; 
(£(E) - £(E+eV))dE,. (2k 4) 


The £(E) can be expressed as 


(Eu) /k,T = 
£-(E) a= ( e + 1 


where E is the total energy for the particular system and 
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yw its chemical potential, which at the low temperatures 
we are dealing with is equal to the Fermi energy. E, is 
the total energy in the direction perpendicular to the 
potential barrier of thickness s and barrier height ¢4(x,V) 
expressed at position x and applied voltage V. 


Brinkman, Dynes, and Rowell (20) entered into this 


expression a trapezoidal model for the barrier 
$(x,V) = $, + (%/s).(b5-eV-9,) 


denoting o> and oo as the barrier heights at metal 1 and 
metal 2 respectively with zero applied voltage. By 
expanding (2.4) in powers of voltage they obtained the 
following result for the conductance o(V). The expres- 
sion, which is accurate to roughly 10% when the barrier 
thickness is greater than 10a and Ao/o is less than one, 


is given by 


A Ad A 
Hale meee O° ets: 2 


= eae = 25 
where Ad = bo-oy i= Zo, +b5), A, 4.(2m)-3 s/3hpiand 
10 -4 roel 
o0(0) =(3.16x10 o,.*./s).exp(-1.025 :346,%).. yExpression 
(2.5) gives a means of obtaining the barrier parameters 
S, b, and Ad (and hence 5 and 5) from experimental data 
provided one knows the junction conductance with respect 


to applied voltage. 
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CHAPTER 3 


INELASTIC ELECTRON TUNNELING THEORY 


Bele LOC LOduUCtION 

Thus far, the tunneling electron has been considered 
to tunnel across the barrier without any energy change 
(elastic tunneling), as indeed it would if the barrier 
were vacuum. But, as the barrier we are dealing with is 
not vacuum but rather a material insulator, there exists 
the possibility that a tunneling electron will interact 
with it in some manner. The questions are: does it, how, 
and to what extent? 

Experimental evidence has shown that the tunneling 
electron does indeed interact with the barrier and alter 
its own energy. This occurrence is particularly evident in 
the derivative of the conductance with respect to voltage, 
(do/dV) versus V, the applied bias, which indicates 
distinctly, changes in the behaviour of the conductance 
due to processes other than elastic tunneling. These pro- 
cesses alter the conductance by only a few percent at most, 
indicating that most tunneling: stiiT occurs elastically. 

The tunneling electron may interact with the 
barrier inmany ways: via phonons, magnons, molecular 
Vibrations, etc. Of interest is the inelastic interaction 
of the tunneling electrons with molecules in the insulating 
barrier. Jaklevic and Lambe (62, 82) first saw evidence 


for this type of interaction when looking at do/dV versus 
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ty 
V spectra for chemically doped and undoped A£-Aluminum 
oxide-Pb tunnel junctions. 

The loss of an amount of energy ev = = Aw, Cora 
vibrational mode with frequency WS of a molecular species 
in the barrier opens up an additional channel for current 
Flow. This results in a break (change of slope) in the 
tunnel junction current-voltage (I-V) characteristic, 
equivalent to a step in conductance versus voltage (o vs V) 
or a peak in do/dV versus V as shown in Fig. 3.1. Because 
molecules can have many different vibrational frequencies, 
and various molecules can exist within the barrier, we 
will see many peaks in the do/dV versus V curve. We call 
this curve the Inelastic Electron Tunneling (IET) spectrum. 
Since the energies (0-500 meV) of the vibrational modes 
we study correspond to infrared frequencies (50-4000 om 1), 
we expect IET spectra to be related to the infrared (IR) 


and/or Raman spectra for the same molecular species. 


3.2. Vibrational Spectroscopy 

Part of the information that we gain with Inelastic 
Electron Tunneling Spectroscopy (IETS) is that of the 
insulating barrier composition. By noting the energies at 
which peaks occur in the spectrum, we can, by the use of 
infrared and Raman data identify the types of molecular 
groups involved. For example, we may see a peak at 360 meV 
a 


in our IET spectrum which corresponds to about 2960 cm 


(wavenumbers). A table of IR identifications will tell us 
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Here, an inelastic process for a vibrational mode 
of frequency Wo results ina break in I vs V, a 


step ino vs V, and a peak in do/dV vs V. 
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that a vibration of this energy belongs to the stretching 
mode of a C-H group. Because IR and Raman identifications 
are important to IETS for the determination of barrier 
composition, these spectroscopies will be discussed. 

First we will consider the nature of molecular 
vibrations and the vibrational frequencies we expect to 
observe. There exist two types of molecular vibrations: 
stretching and bending (or deformation). The stretching 
vibration entails a rhythmical movement along the bond axis 
such that the interatomic distance is increasing or decreas- 
ing. A bending vibration, on the other hand, might consist 
of a change in bond angles between bonds with a common atom, 
or the movement of a group of atoms relative to the remainder 
of the molecule without motion of the atoms in the group 
with respect to one another. As examples, twisting, rocking, 
and torsional vibrations involve changes in bond angles 
with reference to a coordinate system arbitrarily set up 
within the molecules. Fig. 3.2 illustrates various types 
of vibrations. 

In the case of stretching frequencies, assignments 
can be approximated by using Hooke's law. In applying this 
law, two atoms and their connecting bond are treated as a 
simple harmonic oscillator composed of two masses joined 
by a spring. The relationship between the frequency of 
oscillation, atomic masses and the force constant of the 


bond is given by the following equation, derived from 


Hooke's law, 
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motion perpendicular to the plane of the page). 
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where v is the vibrational frequency (oma). c is the 
VeLocity of ligne. (cm/sec), £ is the force constant of the 
bond (dynes/cm) , a4 Me and We are the masses of atoms x 
and y, respectively (g). For single bonds, f is approxi- 
mately 5x10> dynes/cm, while for double and triple bonds 
it is about two and three times this value. It should be 
noted that the frequency is expressed as v = f/c (f has 
Mies sor sec +) and thus has units of cm™+ (called wave- 


numbers) instead of ecu! 


Applying (3.1), to C=H stretching, using SEN ae 


24 


g and i. 64«10 © g as mass values for C and H respectively, 


and 5x10> dynes/cm for the single bond, yields a value of 


: for the frequency of this C-H bond vibration. 


3040 cm 
In actual practice, however, C-H stretching vibrations, 
associated with methyl (-CH,) and methylene (-CH.,-) groups, 
are generally observed in the region between 2850-2960 ae. 
The calculation is not precise since effects arising from 
the environment of the C-H group within a molecule have 
not been accounted for. One often uses the frequency of 
intrared absonption to calculate: the force constants (of 
bonds. 

To use Hooke's law for the calculation of vibra- 


tional frequencies of bond stretching, one must consider 


the relative contributions of bond strengths and atomic 
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masses. For example, a comparison of the C-H group with 
the F-H group, based on atomic masses, leads to the 
conclusion that the stretching frequency for the F-H bond 
occurs at a lower frequency than that of the C-H bond. 
The reverse is true. This is because the increase in the 
ZOLCe COnstant from Lert to. right across the first two 
rows Of the periodic table has a greater effect than the 
mass increase. Hence, the F-H group has a higher vibra- 
tional frequency (4138 cm1) than does the C-H group 
(3040 om). 

One way to study these vibrational modes is by 
passing infrared radiation through a sample of molecules 
and observing the frequencies at which the radiation is 
absorbed. This is called infrared spectroscopy. 

i eae molecuie is to absorb infrared radiation, the 
radiation incident on the molecule must first be of the 
correct frequency to cause a quantum jump in its vibra- 
tional energy. Secondly, the vibrating molecule will 
interact with the oscillating electric field of incident 
electromagnetic radiation if an oscillating dipole moment 
accompanies the vibration. Such a change in dipole moment 
occurs for a molecule whenever a change in position of 
the centers of positive and negative charge resulting from 
atomic motion takes place. When all this happens, energy 
is absorbed and the amplitude of the particular vibration 


is increased. The absorbed energy is subsequently released 
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as heat as the molecule reverts from the excited state 
back to the original ground state. 

One can use a Simple explanation (119) as 
illustrated in Fig.3.3 to show how electromagnetic 
radiation can excite vibrational motion in a molecule. 
First, consider the N, molecule which due to charge 
symmetry has no permanent dipole moment; it can be assumed 
that equal positive and negative charges are present on 
each ae Geneon Now, the electric field of incident 
electromagnetic radiation will cause positive charges to 
move in one direction and negative charges in the other. 
This type of interaction, for a symmetric molecule as No, 
produces no vibrational motion. On the other hand with an 
asymmetric molecule, such as NO, where the negative charge 
on the oxygen atom is greater than that on the nitrogen 
atom, electromagnetic radiation of the proper frequency 
(i.e. the frequency of the electromagnetic wave's oscil- 
latory electric field matches the frequency of the 
particular molecular vibration) can produce a stretching 
of the bond between the atoms as the negative charge on 
the oxygen atom moves in an opposite direction to that of 
the positive charge on the nitrogen atom. This sort of 
interaction results in an enhanced stretching (resonance) 


of the bond between the nitrogen and the oxygen, and 
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ésageqasattys the absorption of infrared radiation. 

Thus, a variation in the dipole moment, u, 
accompanying atomic motion similar to that just described 
for NO allows the vibrational mode to become excited and 
to absorb infrared radiation; such a mode is then said to 
be "infrared active". The intensity of the absorption 
process is found to be greater for increased values of 
the change in dipole moment during a vibration. Quantita- 
tively, the integrated absorption. intensity £Or *a 


particular infrared vibrational band may be expressed as 


SER SO 
307 60 


where N is a constant, c is the velocity of light, u is 
the dipole moment and Q the normal coordinate associated 
with the vibrational mode. For example, the stretching 
vibration of a diatomic molecule has a normal coordinate 
proportional to r-rar where re is the equilibrium inter- 
nuclear distance and r is the displaced value during the 
stretching motion. One can see that the larger the change 
of dipole moment with respect to position during vibra- 
tion, the greater will be the intensity of infrared 
absorption. 

We now turn to the Raman effect which can give 
information about vibrational modes lacking an oscillating 
dipole moment and hence complements infrared absorption 


data. When an electromagnetic wave passes through a 
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transparent material, a certain portion of the radiation 
is scattered by the molecules of the transmitting medium. 
The bulk of the scattered radiation has the same fre- 
quency as the incident radiation, in which case the 
phenomenon is called Rayleigh scattering. The magnitude 
of this scattering is proportional to the fourth power 

of the frequency, therefore radiation at the blue end of 
the visible spectrum is scattered more strongly than that 
at the red end. This fact, taken together with the 
spectral emission curve of the sun, and the human eye's 
wavelength dependent sensitivity, explains the blue color 
of the sky. 

In addition to Rayleigh scattering, molecules in 
the medium can lose energy to or gain energy from the 
penetrating electromagnetic wave. The energy changes 
arise from specific transitions between states of the 
molecules, the excitation of a vibrational mode from the 
ground state to an excited state being an example. This 
is the Raman effect and it occurs as the result of the 
interaction between the optical electric field and the 
polarizability, ofthe molecule. 

To use this technique, one irradiates a sample 
with an intense monochromatic beam (e.g. laser light), 
usually of wavelength 4000-7000 A (visible light). Photons 
of the beam may collide and interact with molecules of 
the sample, and in the process transmit energy to or 


acquire energy from the system (via inelastic collisions). 
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These inelastic collisions result in a modification of 


the photon energy to 


where ue is the incident light frequency, and thy, is the 
energy change in the system where Vv; may be the frequency 
of a vibrational mode. Generally, the total Raman modu- 
lated radiation (at frequencies v57¥4) 1S sabout) 12 of -the 
Rayleigh scattering (occurring at Vo) Also, since the 
vibrational modes of a molecule are more likely to be in 

a ground state rather than an excited state at usual 
temperatures (as determined by the Boltzmann distribution) 
scattering involving energy changes DOV ave) is more 
intense than that for h(v tv). In other words, the photon 
tends to lose energy to a molecule rather than gaining 
energy. in fact for Mee above 700 cmt, scattering corres- 
ponding to the energy change Ley is hardly ever 
observed. 

As an example of how infrared and Raman spectros- 
copy complement each other, consider the co. molecule 
shown in Fig. 3.4. The mode shown in 3.4(a) gives rise 
fOuan Oscillating dipole, and’ will thus be infirared-active 
but not Raman-active. The mode in 3.4(b), however, is 
Raman-active, since a, the polarizability can change 
linearly with the displacement of the oxygens. Since 


3.4(b) does not give rise to an oscillating dipole, it is 
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Illustrated are two vibrational modes of the co. 
molecule. The scissoring (bending) mode in (a) is 


IR-active and the symmetrical stretching mode in 


(b) is Raman-active. 
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not infrared-active. 


3.3. ,Simple Modeis for Inelastic Electron Tunneling 
Having looked at infrared and Raman spectroscopy 
as a means of determining vibrational spectra, we now 
return to inelastic electron tunneling. It has been 
found that when a tunneling electron has sufficient energy 
eV = Aw os it may excite a vibrational mode (with angular 
frequency w) of a molecule in the barrier. Also, since 
the tunneling electron can couple to the dipole moment or 
to the polarizability of the molecule, one may observe 
both infrared-active and Raman-active modes in the IET 
spectrum. The mechanism involved in the interaction is 
Similar to that for infrared and Raman excitations (67). 
In the infrared case one views the dipole moment of a polar 
molecule as having several components - one along each 
bond. These dipoles interact with the transverse electric 
field of the infrared radiation and can absorb a quantum 
of incident sradiation if,the Brequency.is..correct..» In 
inelastic electron tunneling it is the longitudinal 
electric field of the electron which interacts with a 
dipole moment, and it is the electron which gives up the 
necessary energy (E = Aw, Wy is the frequency to be 
excited) to excite the molecule. For the Raman case, the 
polarizability components of molecules are involved. The 
induced dipole moment due to polarizability interacts with 


the transverse electric field of the incident radiation 
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(UV or visible) to produce an excited state. Likewise in 
inelastic electron tunneling the electron's longitudinal 
electric field interacting with an induced dipole moment 
can excite the molecule to a higher vibrational state upon 
absorbing energy from the electron. Thus, the tunneling 
electron can excite both infrared- and Raman-like excita- 
tions. 

The energy level diagram for such an occurrence is 
represented in Fig. 3.5. This figure shows the excitation 
of a vibrational mode by an electron having the threshold 
energy eV = ho.- As the electron energy eV is increased, 
the excitation will still occur but will leave the elec- 
tron some energy. At the temperatures of interest (i.e. 
4.2 K) one can assume that all molecules are in their 
ground states and thus, only processes in which the 
electron gives up energy to excite a vibrational mode 
will occur - the electron will not receive energy from a 
vibrational mode. 

The first explanation for the interaction of a 
tunneling electron with a molecular vibrational mode was 
given by Scalapino and Marcus (110). This was given in 
terms of a simple one-electron picture in which the 
tunneling electron interacted with the dipole field (plus 
its image) of a molecular group in the insulator, near 
to the insulator-metal interface. 

Taking into account the nearest image of the 


dipole, the electron-dipole interaction energy is, using 
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An electron energy level diagram for inelastic 
electron tunneling is shown. Here, a vibrational 
mode of frequency Wo is excited by an electron 


with threshold energy eV = Aw: 
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the dipole approximation, 
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U 3.2) 
where x is the distance from the metal-insulator interface 
normal to the metal electrode to the tunneling electron, 

ry is the component of distance between the electron and 
the dipole perpendicular to the electron's motion (i.e. 

the distance component parallel to the interface), Py is 
the component of the dipole moment normal to the electrode, 
and e is the electron charge. Using a WKB approximation 
for the tunneling matrix element, some other approximations 
followed by the use of the "Golden Rule", and more approxi- 
mations we get a value for the second derivative of current 
with respect to voltage due to the inelastic current as 
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Integrating (3.3) over a voltage increment of a vibrational 
band yields an intensity or conductance change as 
2 
A e f. 


47e m 
> 2m, , a 


pom (3.4) 
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; 2 
where A = sine for a bending mode or A = cos’6@ fora 
stretching mode. Here W, and f., are the frequency and 


oscillator strength of the vth band, 6 is the angle between 


the dipole and the normal to the surface, s is the thickness 
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of the insulating layer, n is the number of dipoles per 


unit area, ¢ =U - Ep (e, 


its value at the Fermi surface and U is the effective 


is the electron energy set to 


insulator potential ignoring the contribution of the 
molecular dipoles), ry LS ae CucoOLL an thervicinity of che 
dipole where the approximation for the potential begins 
to fail - since ry enters into a logarithm, its exact 
Magnitude is not@eritical, and finally (dI/dv) | is a zero 
bias conductance. Considering a hydroxyl group as the 


molecular dipole, and n as 10 OH groups per 100 Re LOL 


oe raking S/X 


the bending mode w, ~ 0.11 eV and f, ~ ro" 
= 30, 6 = 2 eV, and arias = 1/2° the use Of EG... (2.0) Gaves 

a conductance change on the order of 1% consistent with 
experimentally reported results. 

Soon after, Jaklevic and Lambe (63, 82) generalized 
these results to include the interaction of the tunneling 
electron with the polarizability of the molecule. With the 
electron able to induce a dipole moment in the molecular 
species of value p = aE (a varying with molecular vibration), 


and taking into account nearest images, the interaction 


energy is 


re ehe oe (375) 
p (274r 7)? 


Using similar approximation to those of Scalapino and 


Marcus, a value for the inelastic contribution is arrived 
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where N.is.the total number molecules,.B.is a constant, 
and (dI/dv) 9 the zero bias conductance. Integrating to 


obtain the intensity or conductance change one arrives at 


aI 2 
i Pass Og <m|a|0> (3.04) 
(35). L | <m|a|0>| 


with ) |<m|a}o>|? directly proportional to the intensity 
m 
of scattered light as in the Raman effect. Using 


pre em? (from optical Raman data) one gets 


J |<mjalo>|*= 1 
m 
a conductance change of (0.1-0.5)% as compared to 1% for 
the dipole interaction described previously. 

Jaklevic and Lambe also consider the resolution 
of peaks in the IET spectrum, the overall linewidths of 


which are determined by natural widths of molecular 


excitation lines, temperature broadening, and density of 
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states due to the metal electrodes (such as superconductors). 


For normal metals, one may express the peak shape in the 
IET spectrum of a single sharp molecular excitation that 


is thermally broadened as 


V-ii 
dr 7 at (ai=2\e = + (ut2) = = rae 3.8 
me TaRerre)| hat emer. wears Comp iex thecrREAce eee 
av (e--1) 


where T is the temperature and Wo the angular frequency 


of the molecular excitation. This equation describes a 
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symmetrical bell-shaped peak of width 5.4 kT, the height 
of which is proportional to Oana: and the total area 
(intensity) of which is independent of temperature. At 
room temperature, in agreement with experiment, the 
broadening is so great - a peak width of 150 meV for modes 
occurring in a 0-500 meV range - that peak observation is 
practically impossible. For this reason IET spectra are 
taken at low temperatures, such as liquid helium tempera- 
ture of 4.2 K, where the broadening is only about 2 meV. 
The other mentioned factor of importance in 
resolution is the use of superconducting electrodes. In 
this case, one must consider the superconducting density 
of states. Doing this, one finds that the IET peak 
sharpens up, and also for a peak with a sharp natural 
linewidth an "undershoot" is observed on the high energy 
Side of the peak (74, 82, 93). Finally, it should be 
noted that the energy at which a peak occurs is shifted 
from Aw, to (Aw, + 2A) when two identical superconducting 
electrodes are used, where 2A is the temperature dependent 


magnitude of the superconducting energy gap. 


3.4. More Complex Models for Inelastic Electron Tunneling 
After the initial works of Scalapino and Marcus, 
and Jaklevic and Lambe, more complex theoretical treat- 
ments of inelastic tunneling began to appear in the 
literature. Bennett et al. (13), and Duke (39) considered 


inelastic electron tunneling using the transfer Hamiltonian 
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formalism. In this formulation, each metal electrode is 
considered separately, with each described by a distinct 
many-body Hamiltonian, the eigenfunctions of which are 
supposedly known. A third term, known as the transfer 
Hamiltonian, is then added to describe the interaction of 
the tunneling electrons with the molecules or impurities 
composing the barrier. The total Hamiltonian may be 


written as 


with terms for the left- and right-hand electrons as well 
as the "transfer" term. One then views the inelastic 
tunneling process as a transition from a state at one 
metal electrode to another state on the other side via 
interaction with barrier molecules. For the case of 
molecular vibrations, the dipole or polarizability 


potentials previously described may be incorporated into 


the transfer Hamiltonian term, the predicted results being 


in approximate agreement with experimental observations. 
Considering the previous approaches to be 
unnecessarily approximate in nature, Brailsford and Davis 
(19) proposed a stationary-state formulation of the 
problem in which the entire noninteracting one-electron 
system, consisting of an arbitrary interaction potential 
(molecular vibrator, etc.) and a square-barrier potential 


is considered. They consider in detail the special case 
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of a one dimensional square barrier with the vibrator 
interaction within the barrier. In the limit of the first 
Born approximation, the eigenfunctions of the metal- 
insulator-metal system are derived. Since Fermi statistics 
play an important role in the inelastic tunneling process, 
Davis (33) subsequently extended the results to a many- 
electron system. His theory was formally derived for an 
arbitrary tunneling barrier with an arbitrary interaction 
potential between the electrons and a vibrator (ekg? 
molecular impurity, phonon, etc.). The results for the 
inelastic current indicated agreement with the transfer 
Hamiltonian method. In addition, for molecular impurities 
sufficiently close to the barrier boundary, the shape of 
peaks in the IET spectrum are predicted to be fundamentally 
different from those of molecules further inside the 
barrier. For molecules very close to the interface, one 
would see an "anti-peak" or a dip in a*1/av? versus V. 
Schattke and Birkner (16, 112) proposed an alter- 
nate theoretical approach to explain inelastic electron 
tunneling which, like Davis does not employ the transfer 
Hamiltonian approach. Whereas the calculations of Davis 
are restricted to T = OK, and complications can arise in 
the case of many-body interactions for T > OK, the model 
of Schattke and Birkner attempts to avoid these difficulties. 
Their results also predict the same type of lineshapes in 


the IET spectrum for a molecular impurity close to the 


boundary of the barrier, as did Davis. 
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However, they see even more; a distinct asymmetry in peak 
shape and size is predicted for opposite voltage polari- 
ties. Although dips or "anti-peaks" have not been confirmed 
in IET spectroscopy, asymmetry in the "usual" peak sizes 
for opposite polarities has fe Se (75, 44129). 

Other approaches to the tunneling problem have been 
devised@by 1Garokt jet tak... (23), 24,-«25),a26)nn-Combescotiet al. 
(30, 31), and Feuchtwang (46); however, their results are 


not directly applicable to experimental verification. 


3.5. Transfer Hamiltonians since 1976 

In 1976 Kirtley, Scalapino, and Hansma (70) revived 
the transfer Hamiltonian formalism to develop a theory 
for the intensities of vibrational modes in the inelastic 
electron tunneling spectroscopy of organic molecules in 
metal-insulator-metal tunnel junctions. At the outset, 
they bypass the use of a dipole approximation for the 
tunneling electron-molecular interaction, and assume 
instead that the charge distribution within the molecule 
can be broken up into partial charges, with each partial 
charge localized on a particular atom. These partial 
charges originate from the uneven sharing of bonding 
electrons and can be obtained from the dipole derivatives 
of infrared theory. The partial charge analysis allows 
for the description of the interaction at distances 
comparable to interatomic lengths within which the dipole 


approximation breaks down. One can write the total 
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tunneling electron-molecular interaction as 
wie) = eee (3.9) 
3 


where ES and ve are the partial charge and position, 
respectively, of the jth atom. 

The inelastic tunneling matrix element, Megs is 
then calculated using WKB wave functions for the incoming 
and outgoing states, K and K'. The scattering process is 
treated in three dimensions and it can be shown that the 
inclusion of off-axis scattering results in conclusions 
different not only quantitatively but also qualitatively 
from purely one-dimensional calculations. Since the 
electron-molecule interaction is used as a first-order 
perturbation on WKB-approximation wave functions, the 
distortion of the incoming and outgoing wavefunctions due 
to the molecular impurity cannot be taken into account. 
However, the relatively-long-range structure of the inter- 
action decreases the importance of this distortion for 
vibrational-mode couplings. 


The final expression obtained for the second 


derivative of current with respect to electron energy is 


given as 
5 3 27 27 
to, = SES E | ( 5 eye (c' ev) 7 | dg | dg! 
d (ev) A T A fe) fe) 
2 §(fiw-eV) . (3.10) 
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Here, 1 is the current per unit area through the junction, 
n is the surface density per unit area of organic molecules, 
Ep lens) is the Fermi energy in the initial (final) electrode, 
V is the voltage bias of the junction, fiw is the energy of 
the vibrational mode of interest, and L is an arbitrary 
normalization length. Also, it is assumed that the vibra- 
tional mode has a 6-function response in frequency to the 
tunneling electron excitation. The results obtained for 
the peak intensity or change in conductance (obtained by 
numerical integration of Eq. (3.10)) due to inelastic 
tunneling are in good agreement with experimentally 
observed results. One may also predict ratios of peak 
intensities for peaks at opposite bias polarity. 

The theory of Scalapino and Marcus (110) predicted 
that there would be no IETS intensity for a vibrational 
mode which gave rise to an oscillating dipole moment 
parallel to the metal electrode surface. Off-axis scat- 
tering weakens this orientation rule somewhat. A numerical 
calculation based on Eq. (3.10) for the hydroxide ion 
predicts that the IETS peak intensity for the OH stretch 
mode should be 8.8 times stronger when the molecule is 
oriented perpendicular to the metal electrode surface than 
when it is oriented parallel to the metal surface. A 
similar calculation for the infrared-active mode of CO. 
predicts the ratio of intensities for the two orientations 
as 16.4. Thus, although the theory predicts a general 


weakening, the selection rule is still quite strong. 


' does aw 
,epioseionr simapyo Tc sexs. kaw pisces 


,sbouseale (fevi2) Le bhtad aks aa veshas, Lees? ont ae | 

30 yYttens oft ea nuk rs oaBk i RO’. pei gona fev. “ ek 
maersittis ts et J bres saSiesal to Sbois Leno 
-sediy add tsi boeyeds eiial yoriy ) aaprod nity 


yy 
éi4 ot yoneuper? ai sdad@eex nOreaaaTAh S& ee Mascon 
, ~ = 


T62 heniSado etinvea> si .meobterioxs goss <fieisa ‘patde 
Ya -Sentsido) esnssoubsosd net monet: =e vatann cent a 
STesalaqr OF sath (04.6) Ba eT roisexperay fsa. 
VipedaaALase%e is snaiesape boop nt Eis thal 
Aseq Fo eelsbs waives, oats en ond Apa twain 7 ee 
wiitalog said os hBORGE. 48: peered 67 soil 
? puoi acene ste 
i, 
arene ite oa 
a oe 


- 


LStohborg (O%1) SeeorEn HAE 
‘ } 7 2 


noi ehyhaoubwa ead soe rie): aia oe. Bsbed woksal 
tee, Svat ait ted | 
bie ‘waits, ig. 2 aa dasnaadh 
nei sbariue sboxsseup hsm sky 43 dln adregsg 1 
A ,ststave Isdom oft ot Latiarsg bosasixe af at 

gO9 to Sbom svisos+beseatal bead boven ? 
aqoivednehio- ows) esd, ae <ayee oud 4 
Tezaro s etaibedd ¥ | 


detente 80 et? wot yoHen 


6f @ivostom edd ‘nar ? 


4l 


Finally, it should be noted that the above 
theoretical formulation does not take into account the 
polarizability of the molecule. The extra interaction due 
to polarizability may be great enough to be of importance 
and should be calculated. 

The above work calculated the contribution to IETS 
intensities due to a model potential essentially long 
range in nature. A later work by Kirtley and Soven (72) 
peesentsscalctiations offcontributions toyehe*intensitiés 
due to the short range part of the molecular potential. 

In this, the scattering molecular potential is calculated 
using the Xa approximation to the exchange correlation 

potential, and the transfer Hamiltonian approach is used 
to calculate inelastic electron tunneling cross sections. 

It is found that predicted short range cross 
sections are strongly energy dependent and could be quite 
large for tunneling electrons at a molecular bound state 
energy. A comparison, however, of long range and short 
range model predictions for higher harmonic amplitudes and 
opposite bias voltage peak intensity asymmetries indicate 
that the long range interactions dominate in IETS. As an 
example, they consider an A%-AL0,-Rh-CO-Pb tunnel junction 
in which the CO molecule is localized near the Pb electrode. 
The intensity of the CO stretch mode is smaller for AX 
biased positive than for A& biased negative, the ratio 
being ~ 0.68. For A& negative, the electron tunnels 


through the oxide before losing energy to the CO but for 
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AX positive the electron loses energy first and then has 
to tunnel. The asymmetry results since electrons with 
diminished energy are less likely to penetrate the 
barrier Y (7s, Si29), V9 the short range “interaction might —be 
expected to predict a large asymmetry, since the inter- 
action is localized close to one electrode, while a long 
range interaction would give less asymmetry. This is 
just what happens: the current short range model 
calculation described here predicts an asymmetry ratio of 
0.23 while the previous long range theory of Kirtley, 
Scalapino, and Hansma gives 0.59, in better agreement 
with the experimental value of 0.68. The experimental 
observations thus indicate that long range interactions 
dominates inGLeETs. 

Further work using the transfer Hamiltonian method 
for inelastic electron tunneling deals with the orienta- 
tion of adsorbed molecules on a surface. In particular, 
Kirtley and Hall (73) used a dipole-potential in a transfer 
Hamiltonian theory of intensities to study the tunneling 
spectrum of methanesulfonate (CH,SO, ) chemisorbed on 
alumina. Comparison of theoretical predictions with 
experiment supports an orientation of methanesulfonate 
adsorbed on alumina such that the C-S bond is perpendicular 
to the surface. Also, dipole derivatives generated from 
tunneling spectrum intensities compare well with those 
obtained from infrared measurements for similar species. 


In this first attempt to unravel the information 
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contained in all the intensities of the spectrum of a 
Simple molecule, Kirtley and Hall attempted to compare 
3203 : 


They found quite good agreement considering the following 


the dipole derivatives of ethane and propane to CH 


problems: (1) different molecules in different environ- 
ments are compared, (2) polarizability interactions are 
ignored in the theory, and (3) positions of dipoles in 


the molecule may be different from those assumed. 
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CHAPTER 4 


EXPERIMENTAL TECHNIQUES 


4.1. Junction Preparation 

Tunnel junctions were formed by first evaporating 
apcet, Or Circular indium contacts onto a clean glass 
substrate (see Fig. 4.1). Following this, an aluminum 
base film (2000 A thick) was evaporated, and then oxidized 
atiel een: plasma discharge to grow an aluminum 
oxide (alumina) insulating layer. Fabrication was then 
completed by evaporating two desired top metal electrodes 
(2000 A thick). 

The substrate was prepared by cutting a piece of 
pyrex microscope slide 1xl cm square, scrubbing with 
detergent (Alconox), rinsing with copious amounts of hot 
water followed briefly by methanol to quickly remove water. 
Visual inspection of the substrate for freedom from dirt 
and small specks of material followed. A clean substrate 
was then flame polished over a bunsen burner for two to 
three minutes until the corners started to melt anda 
slight yellowish tinge appeared in the glass. It was then 
placed in a substrate holder within the evaporator system. 

Various metals and types of evaporation filaments 
were used as shown in Table 4.1. The metals were loaded 
onto the filaments directly except for the aluminum and 
copper wires. These were first cleaned with a methanol 


rinse and wiped dry with kimwipe tissues. The copper wire 
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Tunnel junction layout. 
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Table 4.1 


Evaporation Materials 


Metal Purity Supplier Evaporation Source 
Ax 59 A.D. Mackay F6-4X.030 W 
cour 
Ag 69 Cominco 539-0055 Ta 
boat 
Au 69 Cominco S39-.005 Mo 
S9A-.005 W 
boat 
Cu electrolytic tough S$39-.005 Ta 
pitch wire boat 
In 69 Cominco S39-.005 Mo 
boat 
Mg 59 Johnson Matthey Si/B=.005° Ta 


covered boat 


Pb 69 Cominco S39-.005 Mo 
boat 
Sn 59 Cominco $39-.005 Mo 


Sao UU Deter 


Evaporation sources supplied by R.D. Mathis Company. 
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(2 mm diameter) was then cut into ~ 4 mm pieces and 
loaded onto a boat filament. The aluminum wire was cut 
into longer pieces which were then bent into the shape of 
Staples and loaded onto a helical tungsten filament. 

Resistive heating of a filament evaporated metal 
onto a substrate into a film shape determined by a mask 
placed right beneath the substrate. The different fila- 
ments, which were mounted on a carousel that could be 
rotated with respect to any particular substrate, had 
cylindrical glass chimneys placed around them to prevent 
mutual cross-contamination of metals during evaporation. 

It should be noted that the substrate holder could be 
rotated if desired and could also be moved upward or 
downward. 

Once the metals were loaded and the desired masks 
were chosen for junction preparation, the evaporator system 
vacuum chamber (enclosed in an 18 inch bell jar) was pumped 
out by a mechanical pump. The lower pressures (~ i602 
torr) used during actual metal evaporation were achieved 
using an oil diffusion pump with a liquid nitrogen cold 
trap backed by the above mentioned mechanical pump. During 
the initial pumping down procedure, the lines of the mani- 
fold allowing for the controlled entry of gases into the 
vacuum chamber, were flushed several times with new gas. 

Some pure oxygen was then let into the bell jar to 
a pressure of 200 microns, and a plasma discharge with a 


current of ~ l1mA and a voltage of ~ 420 V was run for 
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5 minutes. To create the plasma, two circular concentric 
copper electrodes were used, the outer one being a circular 
band of copper sheet and the inner one a circular piece 

of copper rod. The voltage was applied such that the 

outer electrode was positive and the inner one negative. 
This five minute oxygen discharge was the first part of a 
"cleaning" procedure for the system which had been exposed 
to room atmosphere. Oxygen, being reactive is expected 

to combine with stray contaminant molecules which could 
subsequently be pumped out of the system. Following the 
oxygen discharge a longer discharge with argon gas was 
carried out for ~ 15 minutes with ~ 150 microns pressure, 
and 11 mA current with ~ 400 volts. The purpose of this 
was to clean the system by using the relatively large argon 
atoms to dislodge any "dirt" atoms or molecules adhering 

to various parts of the evaporator system within the vacuum 
chamber. The effectiveness of this was reflected Be 
increases in gas pressure (above 150 microns) and in current 
during the course of the plasma discharge. Repeating the 
argon cleaning procedure three times during the fifteen 
minutes allowed for the pumping away of "sputtered" mole- 
cules and the reintroduction of fresh gas. 

Once the cleaning procedures were completed the 
system was pumped down to a pressure of approximately 
ote torr. Indium contacts were then evaporated followed 
by the evaporation of a 2000 A aluminum base layer. The 


film thickness during deposition was monitored using a 
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Sloan DIM=s cuartez crystal monitor. For the aluminum base 
film a deposition rate of 70 Ae was used. A shutter 
system allowed one to establish a desired evaporation rate 
prior to actual deposition onto the glass substrate. 

After the base film was completed the substrate was 
raised to a position 5 cm below the level of the inner 
copper electrode. A plasma discharge technique similar to 
that described by Miles and Smith (98) was used for pre- 
paration of the oxide insulating barrier. The gas mixture 
used consisted of 80 microns of oxygen to which was added 
120 microns of helium for a total pressure of 200 microns 
(403% O5, 60% He). Matheson research-grade oxygen (99.99% 
purity) and helium (99.9999% purity) were used. At full 
output of the power supply (Heathkit PS-3 regulated power 
supply) currents of 2 to 6 mA with voltages of the order 
of 375-430 V were produced. Oxidation times ranged from 
15-30 minutes depending upon the thickness of oxide 
desired. The plasma itself took on a pale ghostly green 
appearance and covered an extensive area whereas a dis- 
charge in pure oxygen was purple pink in color with an 
extent closely confined to the inner copper electrode. As 
shown in Table 4.2 plasma discharge characteristics vary 
as the oxygen-helium ratio is altered. The origin of 
these colors probably lies in the "forbidden" transitions 
responsible for auroral colors due to oxygen (275. 957.2960 


Once the oxidation was complete the system was 


again pumped down to ~ 1x107> torr and top electrodes of 
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Table 4.2 
Plasma Discharge Characteristics of Various 


0./He Gas Mixtures 


Method - obtain desired 0. pressure, then fill to ~ 200 


microns with He. 


Pressure (u) Discharge Color Maximum Maximum 
Voltage Current 
O, O.+He (V) (mA ) 
2 2 
0 200 very diffuse pale blue 410 a O79 
50 230 green ~ 420 me de Sig-0) 
80 210 green 420 3°30 
100 200 ghostly green bluish pink 415 4.0 
oO 200 gentle purplish pink 415 eo 
150 200 purplish pink 415 e635 


185 £35 purple pink 415 ~TL0 30 
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2000 A were evaporated. The first 200 A were evaporated 
Very slowly = in order to disturb the oxide layer as little 
as possible - at a rate of ~ 7 Mee with an increasing 
rate for the remaining 1800 A. 

Once the junctions were completed, the evaporator 
system was vented to room air and the substrate was placed 
in a holder bearing springy, gold coated pressure contacts 
which provided for electrical connection to the junctions. 
Between each springy contact, and evaporated contact of a 
junction, a small flat piece of indium was placed. This 
was done to insure against contact failure in the course 
of annealing studies during which junctions were thermally 
cycled many times and studied for long periods of time (i.e. 
months). After this was done, the electrical ontelintlan Ge 
each metal film was checked; the holder with the junctions 
was then attached to a "dipstick" (see Fig. 4.2) and placed 
in liquid nitrogen. The time elapsed between junction 
completion and immersion in liquid nitrogen was about 10 
Minutes, during which annealing could occur. 

This method of junction fabrication was done without 
exposing the system to the outside atmosphere at any time 
during sample preparation, thus cutting down on contamina- 
ieee Furthermore, once the plasma discharge parameters 
(gas mixture, pressure, discharge current and voltage, 
time, and substrate positioning) were known, a very high 
degree of reproducibility and control of oxidation was had. 


The failure rate of aluminum-aluminum oxide-lead junctions, 
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Dipstick and junction holder. 


f | 


iii 
imu 
\ =o a 
i 


DIPSTICK | 


Dibslick 


DD 


for example, was nil. 


4.2. Modulation Spectroscopy 

Measurements of quantities proportional to the 
junction conductance o (=dI/dV) and its derivative 
do/dv (=d*1/av*) both with respect to voltage (I is the 
junction current) were made using a technique called modu- 
lation spectroscopy. In this, one perturbs a sample with 
a periodic signal (modulation signal) of a certain 
frequency and measures the sample's response (108). 

For example, one may apply a sinusoidally varying 
current (perturbation) of frequency w and amplitude 6 
through a resistor (a linear device). The measured 
response would be a Sinusoidally varying voltage of 
frequency w and amplitude Ré where R is the resistance of 
the resistor. This is shown diagrammatically on the I-V 
characteristic of a resistor in Fig. 4.3(a). Moving one 
step further, consider the characteristic of a device 
shown in Fig. 4.3(b). A voltage response with the fre- 
quency w of the modulation current still occurs but the 
amplitude depends upon the slope of the I-V curve. 

If one now takes a more general curve one may 
obtain the type of responses shown in Fig. 4.4. Approxi- 
mating the above with idealized straight lines one obtains 
signals as shown in Fig. 4.5. These resulting signals 
can be easily Fourier analyzed with the aid of the 


formulas depicted in Fig. 4.6. Hence, the signal in each 
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Fig. 4.3 


Basics of modulation spectroscopy. 


Signal=1I,*5cos wt 
Response = V(I,) +Ricoswt 
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Fig. 4.4 


General characteristic and response to modulation signal. 
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Idealized response. 
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Several signals to be used in the Fourier analysis of the 


idealized response signal are shown. 
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One thus sees that higher harmonics of the fundamental 
frequency exist when one applies a modulated current 
Signal to a device with a non-linear I-V characteristic 
curve. It should be noted that the above Fourier analysis 
merely approximates a real situation, but nonetheless 
gives a mathematical insight into the actual physical 
Situation. 

It is seen that in regions of non-linearity the 
response signal is distorted. Fourier analysis tells us 
that this signal can be represented as a superposition 


of sinusoids of various amplitudes and multiples of the 
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fundamental frequency. A tunnel junction is a non-linear 
device, and so if a current modulation of certain 
frequency is applied, the response voltage will be a 
distorted signal containing frequency and amplitude 
components dependent upon the position along the I-V curve 
of the junction. 

Mathematically, given a modulation current écoswt 
of constant amplitude 6, the voltage V developed across 
the tunnel junction may be written in terms of a Taylor 


series, 


Qu 


2 : 
Vat) ese CL.) qv é6écoswt + LS canals Se Cha cue ae rte 2 
fe) al > 2 + 
e) Ai 


° 
(4.1) 


where V(I) is the DC bias of the junction and w is 27 
times the modulating frequency. One notes that with a 
linear element, (a°v/aI*) | and higher derivatives vanish, 
whereas with the ove wilnase- tunnel junction, they are 
non-zero. Equation (4.1) can be rewritten as 

dca 


v(t) = vir.) +|/ 2% Scoame += Ne ereeonisy 8. 
e) ai 4 aq 
I, I 


(4.2) 


From Eq. (4.2) one sees that if 6 is constant then the 

component of voltage across the qUNnC tion atwangulam.: ne= 

quency w is proportional to (dV/dI), and the component 
O 


2 
ak 20 1s proportional. to (d*v/al Jr Thus, by using a 
fo) 


modulation technique, one is able to produce voltage 
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signals proportional to the desired first and second 


derivatives of the tunnel junction I-V characteristic. 


4.3. Measurement of Tunneling Information 

The dynamical conductance o = dI/dV and its 
derivative, do/dV, with respect to voltage comprise the 
desired information about the tunnel junction. This 
information was obtained by a system using a combination 
Of modulation and bridge techniques with harmonic detec- 
tion, suitably interfaced to a minicomputer facility to 
record and store digital data. 

To understand the principle of harmonic detection 
(1), consider the simplified bridge circuit shown in 
Fig. 4.7. Here, a junction of resistance Ri, serves in 
one arm of an AC Wheatstone bridge while the other arm 
contains the variable resistor Rg set so that Rg x Ry in 
order to balance the bridge; we also set R >> Ra and R,- 
With the bridge thus close to balance, the current through 
each arm is practically the same. With dcoswt being the 


AC (modulation) current in each arm, we have 


Vi = R,Scoswt 


V 


24 


epee 
Pent (dV /dt) Scoswt + 1 (a*V/aI ) §-(ltcos2uc) soe eta 
2 Od Ty 4 I, 


where I, is the DC bias current flowing through the junc- 
S10On,, and Ry ~ Rg at bridge balance. The ditference in 
potential between points 1 and 2, Vi2e has components at 


the modulating frequency w and at higher harmonics. We 
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This figure shows a simple bridge circuit. 


Tunnel 
Junction 


Sweep 


68 


Ow 


aI Sy eri-s not 7 | i iM 


aeq¥ 


‘ 
J 


@ 


7 
Pad 


tOAUW 


see that the signal at the fundamental frequency, Vi2') + 
is proportional to the difference between the dynamic 


resistance (dV/aT) | and R 
fe 


d 


Vi 9 (w) = 6 Beas ae eer coswt. 


The second harmonic component, Vi (2) is proportional to 


ya) : 
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ii 


a pe BAL RES cos2wt. 
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Vi 2 (24) o 
fe) 


These two signals are measured by use of lock-in ampli- 
fiers. 

Direct calibration of dV/di by the use of a decade 
resistance box, enables the computation of the dynamical 
conductance o = dI/dV. In the ideal case, the second 
derivative (a*1/av-) = do/dV is then calculated from the 


second harmonic signal (V5 (20) ) using the identity: 
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The presence of second harmonic distortion in the 
measurement system complicates the above relation and 
necessitates the use of a somewhat different numerical 
method to calibrate do/dV. The details of this are 
described elsewhere (4, 5). We Se able to have 
calibrated values of the conductance o, and its deriva- 
tive do/dV, both with respect to bias voltage. The use 
of a minicomputer based system allows for easy manipula- 
tion of increasing data as well as speedy data analysis. 
Capabilities of the system are described by Adler and 
co-workers (4, 5). Fig. 4.8 shows a current layout of 


the data acquisition and handling system. 


4.4. Annealing Procedure 

The thermal annealing experiments were carried out 
by removing a junction from liquid nitrogen temperature 
and placing it at room temperature, about 20°C, for the 
desired length of time. To prevent the ice which formed 
on the tunnel junction electrodes (when removed from 
liquid nitrogen to room air) from melting and subsequently 
damaging the junction, annealing was carried out ina 
vacuum. After annealing, the junction was again put into 
liquid nitrogen until ready for measurement at 4.2 K in 
liquid helium. 

The voltage annealing experiments were also carried 
out at room temperature in a vacuum. The control junction 


was allowed to undergo thermal annealing (open circuit) 
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while the other junction of the pair was biased at the 
desired voltage and polarity as depicted in Fig. 4.9. 

Since the nature of the annealing study involved 
many measurements on the same junction (i.e. at 0, 2, 4, 
8 hours, etc.) the instrumental settings of the 
measuring system were kept as Similar as possible from 


one measurement to the next. 
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In this set-up, one junction of a pair is voltage 
annealed while the other serves as a control and is 


thermally annealed. 
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CHAPTER 5 


DATA ANALYSIS AND FUNDAMENTAL INFORMATION 


ele LiyCLOGUCELON 

In this chapter, the extraction of basic informa- 
tion from tunneling data will be discussed. Topics include 
peak identification in do/dvV, the IET (inelastic electron 
tunneling) spectrum, determination of barrier parameters 
from Berg ostance curves, and the intensities of IET 


peaks. 


5.2. Peak Assignments 

The mes camene of vibrational modes due to certain 
molecular groups, on the basis of the energy at which 
peaks occur in the IET spectrum will now be discussed. 
Fig. 5.1 shows a typical spectrum of an A&-aluminum 
oxide-Pb tunnel junction. Based on infrared and Raman 
spectroscopic studies, as well as the knowledge of what 
was used in making the barrier, one can identify the peaks 
which occur. Using the work of Geiger, Chandrasekhar, 
and Adler (51), plus that of Bowser and Weinberg (17), one 
gets the following assignments for peak positions in a 
relatively clean A%-oxide-Pb junction as compiled in 
babies. 1. 

The peak which occurs near 450 meV is assigned to 
the O-H stretching of hydroxyl groups located on the 


surface of the aluminum oxide. Since hydrogen was not 
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This figure shows the IET spectrum of an AL-AL oxide-Pb 
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IET Peak Assignments for A&-A’ oxide-Pb Junction 


Peak Position Identification 
(meV) 
37 A& phonon 
Ths) pce OH bend 
90 surface OH bend 
8 AX-O stretch (aluminum oxide) 
Sas) AX-OH bend (of alumina hydrate) 
165 combination and overtone peaks from 75 


and 90 meV fundamental peaks 


239 overtone of 118 meV peak 


450 surface OH stretch 
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deliberately added to the oxygen-helium mixture used in 
producing the oxide barrier, the hydrogen was probably 
due to residual water vapor in the evaporator system. 
Verues Lor OH Stretching via “infrarea data Cild) ®for 
alcohols and phenols are given as 1) 3584-3650 em 2 
(437-445 meV) for the "free" hydroxyl group and 


ave ss00276S0cen — 


(389-432 meV) when intermolecular hydro- 
gen bonding - the interaction via hydrogen bonding between 
OH groups - comes into play. Positions for OH bending 


are given as 1) 1330-1420 em7* (162-173 meV) for OH in- 


plane bending, and 2) 650-769 emt (79-94 meV), a broad 
absorption band shown in the spectra of alcohols and 
phenols determined in the liquid state due to out-of-plane 
bending of the bonded OH group. 

According to Bowser and Weinberg, the absence of 
a peak near 198 meV (1600 om™~) excludes the possibility 
of the OH being associated with molecularly adsorbed 
water (87). They conclude that the hydroxyl groups 
responsible for the 450 meV stretching, and the 75 and 90 
meV bending modes are located on the surface of the 
aluminum oxide insulator (as opposed to being incorporated 
in the bulk as aluminum hydroxide). 

The assignment of the small peak at about 235 meV 
as the overtone of the 118 meV peak is generally accepted, 
however, the identification of the 118 meV peak itself 
still arouses some controversy. Based on the IR work of 


Dorsey (36, 37) with anodic aluminum oxides, Geiger, 
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Chandrasekhar, and Adler (51) attributed this peak as 
being due to the alumina hydrate A%-OH bending mode of an 
aluminum trihydroxide (or trihydrate) even though the 
expected OH stretching mode for alumina hydrate at about 
420 meV was absent. Dorsey suggested that the absence of 
this stretching band could be due to the existence of 
strong hydrogen bonding present within the alumina which 
would reduce the intensity of such a band by broadening 
hts 

Lewis, Mosesman, and Weinberg (85), on the other 
hand, suggests that this peak is due to the A&-O stretch- 
ing mode of an aluminum oxide. The more recent work of 
Bowser and Weinberg (17) strongly supports the idea of 
the 118 meV peak as being due mainly to the Af-O 
stretching mode of aluminum nee, THeEy.-GO, 2n adaqrcion, 
find some evidence for the existence of some aluminum 
hydroxide (monitored by the proposed 118 meV bending 
mode) within the bulk aluminum oxide based on work with 
deuterated samples (as did Geiger et al.). Again, they 
too see no evidence for the AXO-H stretching mode. 

If one looks closely at the region of the surface 
OH stretching peak (Fig. 5.2), one does see a small bump 
at about 415 meV which could be due to the missing A2OH 
stretching peak. The weakness in intensity of this peak 
as well as that contributing to the bending mode would 
suggest that very little aluminum hydroxide exists, and 


that the insulator is composed primarily of aluminum oxide. 
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The surface OH stretching peak at 450 meV is shown 


along with the small peak at ~ 415 meV. 
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It should be noted that the peak at ~ 415 meV could also 
be due to hydrogen bonding of the surface OH groups, or 
due to N-H stretching modes (114). Even so, this does 
not change the conclusion that the barrier is composed 
primarily of aluminum oxide with some aluminum hydroxide 
within the insulator and "free" OH on the aluminum oxide 
surface. 

Other work on aluminum oxide using infrared 
methods reveal some interesting points. Takamura, Kihara- 
Morishita, and Moriyama (120), using anodic aluminum oxide 


ub 


films, observed a peak at 955 cm (116 meV) in the infra- 


red reflectance spectrum of this material. They proposed 
that this was due to an A&-O stretching mode. Slightly 
earlier work of Vedder and Vermilyea (125) observed a 

s (80 meV) 


band at 960 cm. (117 meV) and one at 660 cm 


(transmission spectrum). Both bands were associated with 
AX-O stretching, the 960 em + band corresponding to 
longitudinal modes. 

More recently, Maeland, Rittenhouse, Lahar, and 
Romano (90) studied the infrared reflection and transmission 
spectra of -anodic..aluminum soxide: films)... Thein reflection 
wh 


spectrum showed a band at 960 cm ~ (117 meV) which they 


assigned to the longitudinal mode of the A&-O vibration. 

=] 
The transmission spectrum showed a band at 650 cm (79 meV) 
which was assigned to the transverse mode of the A-O 


vibration. These assignments were based on the work of 


Berreman (14) who suggested that in thin films, one should 
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be-able to observe both longitudinal and OSE Optic 
modes depending upon sample geometry and whether reflec- 
tion or transmission infrared spectroscopy is used. 

These works further support the idea of the 118 mev 
peak as being due mainly to A-O stretching (longitudinal 
mode), and also open up the possibility that one may 
expect to see some contribution from the 79 meV band 
(transverse mode of AX-O stretching). Thus, one has the 
Fo lrowang ¢ 1) tHe -iile meV peak due mainly to A&-O 
longitudinal modes and some alumina hydrate (AX-OH) bend, 
and 2) the broad shoulder by the 118 meV peak due toa 
Superposition of A-O transverse modes of the aluminum 
oxide with bending modes at 75 and 90 meV due to surface 
hydroxyl groups. Resolving the contribution due to each 
mode is no simple problem. 

The last region of interest lies between 160 and 
240 meV, and is made up of relatively low intensity bands. 
These bands are ascribed to overtone and combination 
bands of modes present at 75-118 meV. Bowser and Weinberg 
suggest that the 165 meV peak would correspond to a 
combination of the 75 and 90 meV surface OH bending with 
possible contributions from the first overtones of these 
fundamental modes. As previously mentioned, the 235 meV 


peak is thought to be the first overtone of the 118 meV 


peak. 
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5.3. Data from Tunneling Curves 

Speveral pieces Of information are, taken from the 
tunneling curves. From the conductance curve, o(V) 
versus V, the junction resistance at -30 mV, R(-30 mv) = 
1/0 (-30 mV), is obtained. The region between -30 and 
+30 mV was avoided on account of the large structure due to 
the superconducting lead electrode energy gap and phonons. 

The conductance curve was also fitted with a sixth- 
order polynomial in the energy regions of -500 to -150 
Meveancdero Oto S00 meV. (Fig. 5.3)3 the conductance an this 
region is due primarily to elastic tunneling whereas 
between -150 and 150 meV pronounced inelastic tunneling 
due to aluminum oxide modes occurs. Using the Brinkman, 
Dynes, and Rowell (BDR) equation 7 (20) for elastic 
conductance along with our resistance R(-30 meV) plus the 
linear and quadratic coefficients of the polynomial fit, 
junction barrier parameters were calculated (77). These 
are s: the junction thickness, ¢, the average barrier 
height (¢ = > [o, + o5] at zero bias), and Ago = b5- oy: 
the difference of the barrier heights at the two electrodes. 
Values for the parameters 5 and o> can be easily deter- 
mined from the already calculated ¢ and Ad. The BDR 
model is not taken literally but rather as an indicator 


of the qualitative changes taking place within the 


barrier. 


Finally, the last qualities of interest are the 


areas under peaks of interest in the IET spectrum (do /dv 
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versus V), referred to as peak intensities. From the 
calibrated data, the integrated peak intensity F is 


calculated by 


_ Cd 
F= | & av w) | dav’, (51) 


where On is the conductance used to normalize do/dV (in 
this case oo o(-30 meV)) and g(V) is the smooth back- 
ground under the inelastic peak. The peaks of interest 
were those at 450 meV due to OH stretching as well as 

that at 118 meV due to AX-O stretching (the intensity of 
the shoulder centered at about 85 meV was also found). 

The background for the 450 meV peak was determined in 

the region between -500 to -350 or +350 to +500 (depending 
upon the polarity of the sweep bias which is called 
negative when the A& base electrode is negative) uSing a 
3rd order polynomial as shown in Fig. 5.4. Since the 
choice of a background for the 118 peak was less straight- 
forward, a linear fit using the regions depicted in 

Fig. 5.5 was chosen. The peak intensities were obtained 
by numerically integrating the area under the peak (with 
the background g(V) already subtracted). For the 450 meV 
peak the integration limits are roughly between 385-480 
meV,ewnile: the limits for the 118 and 85 meV peaks are 
shown in Fig. 5.6. A crude estimate of the contribution 


due to each of the 118 and 85 meV peaks is made since they 


are not separately resolved peaks. However, considering 
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the estimates as to background and integration limits 
made for the 118 meV peak, the intensity results changed 
fairly smoothly with annealing and were comparable for 


different junctions, as will be seen later. 
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CHAPTER 6 


RESULTS AND DISCUSSION 


Gare elntroduction 

It is known that tunnel junctions left at room 
temperature for any length of time undergo changes in 
resistance (56). _In order to prevent this, the junctions 
are stored at a low temperature (i.e. liquid nitrogen at 
77 K). To carry out useful studies of phenomena such as 
chemisorption of organic molecules on oxide surfaces, it 
is necessary to have some knowledge of the temporal 
behaviour of the oxide barrier in junctions prior to chemi- 
Sorption. This is “particularly “applicable to cases in 
which one is eee boris chemical changes occurring within the 


barrier (e.g., chemical reactions of molecules infused into 


oxide junctions (64)). 


632.5 Thermal Annealing 

This section describes what happens to A-oxide-Pb 
tunnel junctions during thermal annealing (79). By looking 
at junction resistance, barrier parameters (thickness, 
average barrier height, and the difference between the 
barrier heights atethevAtsand Pb electrodes) andvinelastic 
electron tunneling (IET) peak intensities, one can study 
time-dependent effects occurring in the junction. One 
observes changes in the tunneling barrier and IET spectra 


indicating that surface hydroxyl groups (-OH) experience 
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reorientation. 

The AX-oxide-Pb junctions used were prepared as 
described in Chapter 4.1. One should note that the 
insulating barrier formed in a plasma is not likely to 
be stoichiometric alumina, and may also contain other 
molecular species (such as -OH due to the presence of 
residual water in the bell jar system) besides pure alumi- 
num oxide. The covering metal electrode (Pb) will perturb 
the barrier in at least two ways: the electrode may 
penetrate into the barrier (51, 92) or it may alter the 
Orientation of surface groups, canting them somewhat. Also, 
due to the work functions of the two metal electrodes, a 
strong electric field will arise across the barrier (115). 
For an idea of field strength, a potential difference of 
0.5 V (similar to differences in work functions) across a 
10-A junction (these are typical barrier thicknesses found 
using ellipsometry (80)) gives a field of 5x10° V/il. Unis 
field can aid in group reorientation and ionic diffusion. 

The thermal annealing experiments were carried out 
by removing a junction from liquid nitrogen temperature and 
keeping it at room temperature, about 20°C, for the desired 
length of time. To prevent the ice which formed on the 
tunnel junction electrodes (when removed from liquid 
nitrogen to room air) from melting and subsequently damaging 
the junction, annealing was carried out in a vacuum. All 
annealing times referred to are cumulative times (i.e., the 


total time for which the junction has been annealed). 
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The most obvious parameter which changes with 
annealing is the junction resistance R. Fig. 6.1 shows a 
plot of the resistance ratio RL/R, versus t (total 
annealing time) for a series of thermally annealed (open 
GirCuit)) Oxide junctions where Ro is the initial resistance 
(within a few minutes of junction formation), and Ry the 
besistance at time t (at =-30 mV). The curves suggest three 
regions: (1) an initial decrease (some curves drop below 
Ri/R, ce O0)melneresestance (or none at all) >. (2). alrapid 
resistance increase which becomes (3) a more gentle 
increase. 

Other noticeable changes occur in the conductance 
‘oO versus applied bias V, and in the IET spectrum do/dv 
versus V. The spectrum of an oxide (Fig. 6.2) shows two 
main peaks of interest, the 450 meV peak due to surface OH 
stretching, and an AL-O stretching mode at 118 meV due to 
the aluminum oxide (all peak positions referred to are for 
the polarities AX negative and Pb positive). Upon thermal 
annealing one sees noticeable changes in the spectrum of a 
Partaculan jJunetion: “for example (Fig, 623) ,) the OH peak 
at 450 meV increases in size and sharpens up slightly. 

Changes also occur in the shape of the conductance 
curves. These curves are characterized by the barrier 
parameters (see Chapter 5.3) s, the average barrier thick- 
ness, $, the average barrier height [¢ = 5 ($5 +65) at zero 
bias], and Ad = d.-¢,, the difference of the barrier 


heights at the two electrodes. For thermally annealed 
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Feige. 
This figure shows a plot of R,/R, vs t for several 


junctions, where R, is the junction resistance at 


c 


time t, Ro is the initial resistance and t is the 


total annealing time. 
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A typical IET spectrum for an A&-Af oxide-Pb tunnel 
junction is shown. Peaks of interest occur at 118 


meV and 450 meV. 
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The change in the 450 meV OH peak with annealing is 
shown here for two different times - 0 hour and 
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Junctions, one finds that s, the effective barrier 
thickness, increases with time as shown in Fig. 6.4 (this 
does not mean the barrier is necessarily becoming physi- 
cally wider but rather the tunneling resistivity is 
becoming greater), and that $, the average barrier height, 
G@ecreases with time. The most interesting results involve 
the barrier heights at the A and Pb electrodes, 5 and 5 
respectively. Figure 6.5(a) shows a decrease in Ad which 


teveils Ort within 40 hours. In Pig. 6.5(b) one sees that 


Oi" the barrier height at the aluminum electrode, decreases 


and then levels off. This suggests that there has been 
an increase in positive charge density near the lead 
electrode. In conjunction with this barrier height 
behaviour, one finds that the intensity F(450) of the 
450 meV surface OH stretching peak for this junction 
increases with time and tends to saturate. 

Looking back over the results we find 


(1) a rapid 450 meV peak intensity increase to 


saturation, 
(2) a decrease in Ado, 
(3) a decrease in b> (the lead side of the 


junction becoming more positive), 


(4) 1 remains fairly constant for the aluminum 


electrode, 


(5) a slight sharpening of the 450 meV peak. 
Since a large electric field may exist across the 


barrier (~ 10° v/m) a possible interpretation can be 


106 


~gOT 


ei yt ivedutesy cat Seine al eons: tod e 
ttoiet ssivtad sesvave ats e ‘any tea: nes 
eviovar exiuv2s5 nia kemelhias Shadlbnaid purets ied ' 
has -?) yesbermootn ak fin aK oe da 4 pee ne | 
foldw A wt eaessceh s eoade- (s7ere omit a 5 oh . 
Jad soeq. ono (a}2.8 ode a eerie! 7 
zezagineb ,sborisele munities a a unmae 
tel suede sede ereogune east s Xho: te w TS, 

set $03 sen yitersh open eet laecit alia a 

weaien s9itzsd ehele rh the nékiPematies a2 dt, 

ai? 4 ine Ny {2 SB ria ea ts 

: sigh 208 ask Ata RO -ag 
oFeusie os (Sti ‘oe nat Atbw i 

Gus? ew Sohbet ona pone sed pakecct. 
at seseronti. vodeseans Riess vem Geb REgs* a eu 


» 


Nels Say 


a) eK seeerseb se ft) 
aa: Fo. shed Saal qa) ae Wi sasesseb s (8) 
; «(orks Ealag ‘piom pabmosed’ ol senit 
sunimeis 21% 103 as tre ‘babes watemex b () 


Z etoutpele : : 


Riser Vert 022 site 38 aiiaid sess m  (@) 


silt aanios taine yen bled QEeemele epued # oenle 
od as? nodes sidyecet) aldivehg s takv Cos -} <eivisd 


=) 
é 


\ 
jc 


ere int hua oii Si © sisiw ,mol 
shea te da “lieve! (ie aes SF Son SLuoudn ej 
skBeed en wakyaluptes xine x04 % ¢¥ed ae te} 
im | . -te.5)) pe 


Fig, son4 
A plot of s vs t is shown for a thermally annealed 
junction, where s is the barrier thickness. These values 
of s should not be taken too literally as they are 


influenced by ¢ for their calculation as indicated in 
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Changes in Ad, da: do: and F(450) are shown. (a) ) A 
plot of the separation between barrier heights at 
the two electrodes, Ad vs total annealing time t 
shows the decrease of Ad with time. (b) This plot 
of the barrier heights at the two electrodes vs 
time shows a decrease at the Pb side ($5) but 
little change at the A’ side (o,)- (ole eihias “aap 
shows how the 450 meV peak intensity, F(450) increases 


with annealing time and finally levels off. 
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given by considering the reorientation of the surface 
hydroxyl groups on the aluminum oxide. The above 
mechanism is shown schematically in Fig. 6.6(a) where the 
originally quite randomly aligned OH groups reorient 
themselves under the influence of the electric field. 
Fig. 6.6(b) shows that, in realigning the OH dipole, 6 
(the average angle of an OH group from the normal to the 
plane of the Pb electrode) tends to 0°, i.e. the OH 
stretch motion is more perpendicular to the Pb electrode. 
This will result in an intensity increase for this mode 
proportional to cos“6 for a stretching =mode (110). 
Secondly, Fig. 6.6(b) shows that with reorientation the 
positive hydrogen of the OH moves closer to the lead 
electrode, increasing the positive charge density near it; 
therefore 5 decreases while practically no effect is 
exerted on 5 at the A& side. 

One also finds that upon annealing, the intensity 
of the 118 meV peak due to AX-O stretching increases in a 
manner similar to the 450 meV peak increase. By looking 
at the intensity-vs-annealing time curves (Fig. 6.7(a) for 
a few junctions, one can see that the lower the initial 
118° mev intensity the greater the annealing effect. 
Whereas the 118 meV peak intensities tend to saturate at 
different values, the 450 meV intensities all tend to 
saturate at approximately the same value feid.26. 0 bie 
This suggests that there is a certain "monolayer" coverage 


of surface OH for all junctions (each of which has its own 
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Rig .4.6%50 
(a) A schematic representation of the tunnel barrier 
shows how the OH groups reorient under the influence 
of an electric field during thermal annealing. 
(b) Here, a close-up view shows how 6 decreases with 
reorientation and results in the hydrogen side of 


the OH group moving closer to the Pb electrode. 
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(a) A plot of the 118 meV (A£-O) intensity vs 
annealing time for several tunnel junctions. 
(b) This graph of the 450 meV (O-H) intensity vs 
annealing time for several tunnel junctions shows 
that for all the samples the same saturation 


intensity is reached. 
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particular initial orientation) while the amount of 
aluminum oxide can vary [Fig. 6.7(a)]. 

To summarize, it is found that A-oxide-Pb 
tunnel junctions subjected to thermal annealing undergo 
resistance changes (up and down) as well as changes in 
barrier parameters and IET peak intensities with time. 
This suggests a reorientation of the surface hydroxyl groups 
under the influence of an inherent electric field of polar- 
Poye Such? that the Atssidelas positive ‘and the® Pbh'side?is 
effectively negative. Also, by looking at the values at 
which the IET peak intensities saturate, one deduces that 
there is a certain monolayer coverage of surface OH for all 


the junctions while the amount of aluminum oxide varies. 


ba3cee Voltage’ Annealing 

Chapter 6.2 has described some of the effects which 
occur in Af-oxide-Pb tunnel junctions when they are 
thermally annealed (at room temperature). This section 
Gescribes\ the. restilts#’of voltage annealing “(80).*)"™inSthis, 
arvconstant® voltage is’ applied) across the-jgunction elec 
trodes to produce an electric field which adds to or sub- 
tracts fromuthe: junction’ s intrinsic electrucy freldee This 
intrinsic field is responsible for some of the changes 
observed in thermal annealing; thus, altering the field 
with an applied bias should help to clarify our knowledge 
of the annealing process. Information about junction 


resistance, barrier parameters, and inelastic electron 
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tunneling (IET) peak intensities indicates that positive 
charge at the A&-oxide interface can be increased or 
decreased, depending upon the polarity of the applied bias. 
Also, by proper choice of bias polarity, junction aging 
(i.e. any changes occurring within the junction after its 
completion) can be retarded. 

The A£-oxide-Pb junctions used were prepared as 
previously described such that two "identical" junctions 
Were made.On each substrate. .;in this.way ,,onesjunction of 
the pair served as a control (thermally annealed only) 
while the other junction was subjected to voltage annealing. 
The voltage annealing experiments were carried out at room 
temperature in a vacuum. The control junction was allowed 
to undergo thermal annealing (open circuit) while the other 
junction of the pair was biased at the desired voltage and 
polarity. The biasing arrangement with the A& electrode 
positive, Pb negative shall be referred to as forward 
polarity. (FP) and AX negative, Pb positive as reverse 
polarity (RP). Previous thermal annealing experiments 
showed that the intrinsic field of a junction was of forward 
polarity. 

Again, interpretation of the annealing results is 
carried out by examining changes in junction resistance, 
barrier parameters, and inelastic electron tunneling peak 
intensities with respect to total annealing time. In 
addition to those exhibited by thermally annealed junctions, 


voltage annealed junctions showed several new effects. 
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Junctions biased at 0.5 V with forward polarity 
(AX positive, Pb negative) show an increase in the slope 
of the resistance versus annealing time curves shortly 
after the desired voltage is applied. In Fig. 6.8 a 


graph of R,/R, versus time t-t, (where R, is the junction 


e 
resistance at time t, Ry is the resistance at the time of 
bias application the and t is the total annealing time) 
illustrates this resistance increase for several junctions 
of different ages when initially biased. One sees that 

the increase soon slows down again and that the net resis- 
tance ratio increase after, say 50 hours of application 

of the bias, diminishes the older the junction was to 

start with. 

Another noticeable change 1s that which occurs’ in 
oy and Por the barrier heights at the AX and Pb electrodes, 
respectively (depending on the nature of the sample, 5 
can be greater or smaller than d5)- Following application 
of an annealing bias, it is found that Ad = d5-o) increases 
in time and finally levels off as shown in Fig. 6.9 by 
Ag-Ad, versus t-t, (Ady refers to the value of Ad at 
time th): It is seen that an age dependence is also 
exhibited by the magnitude of the effect. One should note 
that: it is ory which decreases markedly while 5 changes 
only a little. This implies that the insulating barrier 
at the aluminum electrode acquires positive charge with 
time which eventually saturates as long as the biasing 


voltage is maintained. This charge, moreover is trapped 
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This figure shows a plot of RL/R, vs t-t, for three 


junctions biased at 0.5 V FP. RL is the junction 


resistance at time t, Ry is the resistance at the 


time of bias application t ana t-is -the total 


lend 
annealing time. For the circles, squares, and 


triangles ty = 0, 212, and 324 hr. respectively. 
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ADO tarot Ag-Ad, versus t-th for three junctions 
shows the effect of a 0.5 V FP bias on Ad, the 
difference of the barrier heights at the two 
electrodes. The junction age when biased, th. 
equals 0, 212, and 324 hr. for circles, squares, 
and triangles respectively. Ady is the value of 


Ad at time th: 
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just within the barrier otherwise, like a Gapacitor, Lt 
would drain off as the annealing voltage was reduced to 
zero. 

If junctions are annealed using reverse polarity 
(AX negative, Pb positive) at 0.5 V, the results are 
aitferent. Looking at R,/R, yereueae (Ro is the resistance 
at t = 0), one finds an. initial decrease in resistance 
which eventually increases slowly (Fig. 6.10). In addition, 


the barrier height, $6 at the A& electrode goes up 


1’ 
(implying that the AX side has become less positive) while 
that at the Pb electrode, Por remains constant. 

These resistance changes as well as the behaviour 
of 1 Crg036..4 1) and 5 for a biased junction suggest 
movement of positive AX ions into and out of the barrier 
at the Af-oxide interface. When forward polarity is 
utilized, aluminum ions may be created at the aluminum 
electrode and drift a short way into the oxide barrier 
under the influence of the applied field. Such an occurrence 
raises the amount of positive charge near the AX electrode 
(do) thus decreases) and increases junction resistance since 
the inclusion of AX ions within the oxide makes the barrier 
denser near the A&-oxide interface. On the other hand, 
application of a reverse polarity bias causes an Lnwtaal 
resistance decrease and an increase in di: This indicates 
that some already existing positive charge is moved out of 


the barrier, which is reasonable since a negative potential 


has now been applied at the A% side. 
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Figi, s017L0 
The resistance behaviour, RL/R, versus t is shown 
here for a junction biased atv0, owe el (cineles):, 


and a thermally annealed control junction (Squares). 
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This figure illustrates the behaviour of the barrier 
height at the AX electrode, 4 vs time for junctions 
subjected to 0.5 V. RP (cirelles)770 (5 V-Fe (squares), 


or thermal annealing only (triangles). 
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One can account for the age dependence of the 
amount of change in R,/Rp» and in Ad by noting the 
behaviour of a voltage annealed (FP) junction once the 
annealing bias is removed. Fig. 6.12 shows what happens. 
A decrease in Ad [Fig. 6.12(a)] occurs and o1 increases 
Suggesting that the amount of positive charge near the AX 
electrode is being depleted. Looking now at the junction's 
resistance curve [Fig. 6.12(b)], it is seen that as charge 
depletion occurs, resistance increases. This leads one 
to believe that AX ions are being neutralized, and remain 
within the barrier as opposed to drifting back out to the 
AX electrode. The age effect can now be explained as 
follows. As a junction anneals, A ions are created and 
deritwinto che Junction “due ttowles intrinsic. field or under 
the influence of an externally applied field. The ions, 
once having entered the barrier are eventually neutralized 
and new ions can be formed (maintaining a type of dynamic 
charge equilibrium). As a junction ages, vacancies for 
new ions to occupy within the insulator are diminished since 
some have already been filled. Thus, when an external bias 
is applied to an old junction as opposed to a newer one, 
its response is not as great simply because there are not 
as many vacancies in the barrier for new A& ions to fill 
and produce their effects. This is illustrated in Fig. 6.13. 
On the subject of aging, one finds that junction aging can 


be suppressed by the application of a reverse polarity bias. 


This will be discussed later. 
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(a) Here we see the decrease in Ad after the 
removal of a 0.5 V FP bias. The arrow indicates 
the beginning of thermal annealing. 
(b) The resistance ratio behaviour, R,/Roy Ores 
Aeryerenon during voltage annealing (0.5 V FP) and 
then thermal annealing (starting at the arrow) 


is shown. 
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(a) This figure shows the response of a “young” 
junction to FP biasing. In the first picture many 
vacancies exist at the AX-AX oxide interface. When 
a bias is applied (second picture) ions are created 
and can drift into the barrier to fill the vacancies. 
They remain trapped here after the bias is removed 
(picture three). (b) Here, the response of an 
"old" junction to FP biasing tis=shown.. | Vacancies 
are filled in the same way as described in (a) 
except that there are fewer vacancies to begin with. 
Circles with + signs represent A& ions, and shaded- 


in circles denote A atoms. 
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Now, one turns to see what effects voltage 
annealing has on IET peak intensities. Looking at an 
aluminum oxide spectrum (Fig. 6.14) one notes the usual 
peaks of interest at ~ 118 and 450 meV due to A2-O and 
O-H stretching respectively. When voltage annealing is 
Carriea Out with a 0.5 V FP bias, one finds the 118 and 
450 meV peak intensities increase with time to values 
greater than those in the thermal annealing case. This 
suggests that with an electric field greater than the 
intrinsic field alone, the surface hydroxyl and aluminum 
Oxide may experience further reorientation. On the other 
hana, thie application of a field of 0.5: V RP is found £0 
inhibit the change in peak intensities which normally 
occurs in virgin junctions under the influence of their 
Ppoteinsic £ileld (Figs 6.15). > this amplies that the 
intrinsic field of the junction is neutralized by the 
0.5 V RP bias, thereby leaving no electric field across 
the junction to reorient the molecular groups within. If 
the bias is now removed, the junction thermally anneals 
as before; resistance and peak intensities increase. By 
applying such a reverse field to neutralize the intrinsic 
field, one can thus delay the aging of a tunnel junction. 

Having looked at the voltage annealing behaviour, 
one can now use it to aid in explaining some results 
observed in the thermal annealing experiments. A number 


of thermally annealed junctions initially showed an 


unexplained drop in junction resistance Similar to the one 
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Pigz: 1o.el'4 
The two curves are IET spectra for an AL-AX oxide-Pb 
junction at two different times. The upper curve 
Shows the effect that 70 hours of 0.5 V FP voltage 


annealing has on the virgin sample (lower curve). 


- 4918118 H-0 OSp 


ie) 
QO_ 
' 
w 
BS 
x 
Oo 
§ 
< 


Y9}911S Q-|V 8ll- 


E (meV) 


1S5 


a a 


E (46 Ay 


ae. Ly i a : 
7 ; : 7 af 
Were : 
fe 7 . 
oe ae Ty : fi A 
7 Hb cs / 7 - ae : a 136 “4 
LP Jip 7 


b emis egatov 4 swatgansak freq Fo} aly eiaT 


7 rei Maomedtqqne (Vo 2, OF» vate vs iwalee Fader: 


ae edd ME. exwon OSL = yp movd) to sepneto 
a eiseeds BoA ate «Af ent 4 bar OUTS s 
soupigront ers peau 
po aa, Sates 
i 
ee | a 
} ’ 
ee ip « -- 
Z 
‘ “ " ' 
Th , 2 
8 o is 
* 4 : s 7 
1 


EA Gi Omics 
This plot of peak intensity, F versus time shows how 
reverse polarity annealing (0.5 V) suppresses any 
changes in F (from 0 - 120 hours). If the bias is 
removed (at the arrow), the junction thermally 


anneals and peak intensities increase. 
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observed in voltage annealing with reverse polarity. 
Unlike RP voltage annealing the peak intensities increased 
in the usual way. A close look now at 5 for those 
junctions with a pronounced resistance dip shows that 1 
initially increases a little (A side becoming more 
negative). This, as well as the resistance drop, could be 
eecounted fonei ft positive AY ions lett the’ barrier. 
However, Since the intrinsic field of the junction has the 
AX electrode positive, there seems to be a contradiction 
unless one considers what may occur in the course of 
Bunecion febrication. Looking at. the barrier during 
growth, one might imagine it to be aluminum oxide with 
positive AX ions near the aluminum electrode and OH groups 
standing up on the growing alumina surface [Fig. 6.16(a) ]. 
Subsequent evaporation of the Pb cover electrode will tend 
to squash the OH groups, leaving them canted as shown in 
Fig. 6.16(b). This movement of the OH group's hydrogen 
closer to the oxide will cause a greater electrostatic 
repulsion between the aluminum ions and the hydroxyl 
hydrogens than existed prior to the deposition of the Pb 
cover electrode. Since the A& ions are smaller and more 
mobile than the OH groups, the mutual electrostatic 
repulsion can force some ions out of the barrier 
[Fig. 6.16(c)]. This takes place only until the intrin- 
sic field of the junction can reorient the OH groups 
sufficiently to relieve the repulsive force. At this 


point, A& ions under the influence of the intrinsic field 
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Pig 26-16 
This figure shows how evaporation of the Pb electrode 
on the growing oxide can perturb the surface 
hydroxyl groups and cause expulsion of A ions out 


of the barrier due to electrostatic repulsion. 
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are once more able to enter the barrier [Fig. 6.16(d)] 
and cause resistance increases. Some thermally annealed 
junctions show the above behaviour; some do not. Factors 
influencing this probably include (1) the manner in which 
the Pb electrode is evaporated and (2) the number and 
grcstribucionvof AY ions within the oxide, 

In this section, a new technique for the study 
of ion motion and behaviour within tunnel junctions has 
been introduced. By applying a 0.5 V FP annealing 
voltage to a junction, one causes positive charges, pre- 
sumably aluminum ions to accumulate at the A&-oxide 
interface. These ions serve to increase Ad (do, decreases) 
as well as the junction resistance. When the voltage is 
removed, neutralization of ions within the barrier is 
reflected by the fall of Ad Co) now increases). The 
accumulated charge also serves to enhance the junction's 
intrinsic electric field which further aids group reorienta- 
tion. On the other hand, when a 0.5 V RP bias is applied, 
some positive charge is removed from the barrier and 
junction!lresistance initially @allste Tternises slowly 
afterward, perhaps due to neutralization of some remaining 
charge within the barrier. Peak intensities do not change 
because the magnitude of the bias seems sufficient to 
negate the junction's intrinsic field, thereby leaving no 
electric field to reorient molecular groups. Junction 
aging is thus postponed until the reverse polarity bias 


is removed, at which point normal thermal annealing occurs. 
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Finally, an explanation for the resistance dip Hee or 
in some thermally annealed junctions has been suggested 
as due to a temporary electrostatic repulsion occurring 
between A& ions and the H of surface hydroxyl groups 


after the Pb cover electrode has been evaporated. 


6.4. Voltage Annealing with Alternating Polarity 

Having studied annealing with an applied voltage 
of given polarity, attention is now turned to the effects 
of voltage annealing in which the polarity of the applied 
bias is reversed for each annealing period of time ey 
Such annealing is referred to as voltage annealing with 
alternating polarity. This is done in order to see 
whether any of the annealing effects brought about by a 
voltage of given polarity can be reversed by applying a 
voltage of opposite polarity for the same duration of 
time. Although more work needs to be done, some results 
can be shown. 

For this experiment we use two Af-aluminum oxide- 
Pb tunnel junctions grown on the same aluminum base layer. 
The tunneling characteristics of the virgin junctions 
are firstemeasured.  sFollowingsgthis, theycontroljis 
thermally annealed, while the other junction is subjected 
tosatbias.of either polarity for a time As in the 
previous section, we refer to the polarity with the A% 
electrode positive, Pb negative as forward polarity (FP) 


and AX negative, Pb positive as reverse polarity (RP). 
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At the end of an annealing period, the two junctions are 
cooled and their tunneling characteristics are again 
measured at 4.2 K in liquid helium. Changes occurring 
in the tunnel junctions are studied with respect to the 
total annealing time, t. 

Previous results with voltage annealing (80) showed 
that the application of a room temperature annealing 
voltage of forward polarity could force positive charge 
into the junction at the Af-AX oxide interface, some of 
which was neutralized within the junction. Also, a 
reverse polarity voltage was able to draw a small amount 
Of positive charge out of the barrier. This, along with 
other information suggested that some of the charge was 
fairly mobile - it could be moved in and out of the 
barrier easily - while another portion became trapped 
eee the Cranes barrier. 

To explore this further it was decided to apply a 
voltage of forward polarity to a virgin junction for time 
t_ in order to move charge at the aluminum electrode in 
toward the junction interior. After that time a reverse 
polarity voltage would be applied for time ce again to 
try and counteract the effect of the previous voltage 
anneal, thereby moving charge back to the aluminum 
electrode. 

Figure 6.17 shows the resistance behaviour of a 
junction for which this has been done. The voltage 


magnitude was 1 volt and the annealing period t was 2 
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FiGen6..M, 
The resistance behaviour for an alternating polarity 
voltage annealed junction. Here, the magnitude of the 
voltage is 1 volt and the annealing period eS isi 2 hours: 


The thermally annealed control is shown below (squares). 
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hours. One sees a large resistance increase in the first 
two hours. Application of the reverse polarity voltage 
causes a comparatively small resistance decrease, while 
subsequent application of the forward polarity voltage 
again, causes a small resistance increase. As can be 
seen, resistance goes up or down depending upon the 
polarity of the applied annealing voltage with a general 
trend of overall resistance increase. One also notes that 
the size of the up or down swing decreases with time. The 
thermally annealed control is also shown. 

In Fig. 6.18, the behaviour of s, the calculated 
barrier thickness is shown. The first 2 hours exhibit a 
large change in thickness while subsequent changes are 
smaller, and again increase or decrease depending upon 
the annealing voltage polarity. 

Fig. 6.19 illustrates the behaviour of 5 and d5- 
Looking first at dys the barrier height at the Af-aluminum 
oxide surface, we see that with a forward polarity bias, 
5 drops ‘dramatically in the fist two hours indicating 
an increase of positive charge at this interface. Subse- 
quent application of a reverse polarity bias causes or 
to increase somewhat, suggesting a decrease in the amount 
Of positive charngerat the interface. — further annealing 
with opposite polarities again results in an oscillatory 
behaviour in o4 with a magnitude of the up or down swing 
decreasing with time. Similarly, bo also exhibits an 


oscillatory behaviour. If one includes the results for 
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Pilg. Oc 
The change in barrier thickness, s, with respect to 
total annealing time, t, is shown. Again the 
magnitude of the voltage is 1 volt and . is= 2 hours. 


Squares represent the control junction. 
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Fig. siOr Ho 
The behaviour of aT and 5 (barrier heights at the A’ 
and Pb electrodes) versus time for alternating polarity 


voltage annealing is shown. 
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ib saih 
Pnescontrol junction (Pig. 6.20) 2t as seen: that ‘the 


control values of b> - the barrier height at the lead 
electrode - lie within the region determined by the 
oscillating b> values of the other junction. We see 
though that the same is not true for oi the alternating 
polarity has resulted in a large downward shift in os: 
This implies that the interesting things are happening at 
the AX electrode. 

Summarizing this data, we see that there is a 
large change, occurring at the AX-A& oxide interface, in 
the first 2 hours when the junction is subjected to an FP 
anneal which cannot be completely reversed by an RP 
anneal. Moreover, subsequent changes due to FP annealing 
are small compared to those in the first 2 hours. This 
suggests that initially, the influence of the FP voltage 
quickly forces a portion of the positive charge, presumably 
AX ions, into a number of available trapping sites from 
which the ions cannot be easily removed when a reverse 
polarity bias is applied. By trapping is meant any mechan- 
ism within the barrier which makes it difficult to later 
remove the charge from the barrier. Such trapping could 
occur by the filling of vacancies, or by bonding of Av 
ions to negative ions within the insulator, possibly oxide 
ions. 

The results following the first two hours show 
that a smaller portion of charge seems to be fairly mobile 


under the influence of changing polarities, however the 
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Hig2aOraa0 
The behaviour of oy and 5 for the thermally annealed 
control junction is included in this picture as denoted 


by the shaded in circles and squares. 
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amount of this mobile charge seems to be decreasing with 
time as indicated by the decreasing amplitude of up and 
down swing on the resistance and oy graphs.-<# This along 
with the general increase in resistance suggests that 
chaege. trappingsis.stilld occurring but ate a.much-slower 
mate.than in.the initial.2 hours. 

1O,0btain?y more;ansight, into, this,work,efurther 
study should examine the influence of such variables as 
1) length of annealing period, tf 
2) magnitude of the annealing voltage, and 
3) the effect of initially annealing the junction with a 
reverse polarity instead of forward polarity. 

One should note that in the above experiment the 
intrinsic. fieid ofsthe.junctaoneiseanot consideredxenThis 
field has been shown to be of forward polarity with a 
magnitude,of,about,s0..5 volts. Thus, “applying apbias .of 
1V FP yields a net field of 1.5 V FP while a bias \ of 
LVoRP.givesija.net.field of 0.5 V RP~ Therefore, if one 
wishes to have equal and opposite net fields across the 
jJunetion of 1 volt magnitude, one must use asymmetric 
applaedmbiases @f O0254VeEP and 1.58V Reaeespectively= york 


has.begun-on this but, the results are#mot yet »complete: 


6)544O0nigin Ofathesintrinsic Electhexcehienhd 
At this stage, the origin of the»junction's»intrin- 
Sic electric field is considered. The intrinsic field of 


the AL-A& oxide-Pb tunnel junction probably does not arise 
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solely from a difference in the work functions of aluminum 
and lead (AX = 4.28 eV, Pb ~ 4.25 eV) since this difference 
mo. qutte small )(97)) compared to the estimated 0.5° V 
potential difference. Also, a field seems to exist in 
AX-AX oxide-A& junctions for which the two electrodes are 
the same, thus it would seem that something within the 
oxide insulator is primarily responsible for the intrinsic 
electric field. This “something” could be the presence 
of a relatively immobile (due to the ion size) negative 
oxide ion layer created in the process of oxidizing the 
aluminum. It would be near the surface of the oxide (i.e. 
closer to the Pb electrode than the A one in an A-AgX 
oxide-Pb junction) and thus create an electric field of 
forward polarity (Pb-ve) as observed. The presence of such 
oxide ions has been described by Mott (102) in his theory 
of the formation of protective oxide films on metals. 

Such a negative oxide ion layer (Fig. 6.21) could 
have several effects: 
1) it would provide a force for the movement of positive 
AX ions in the barrier; 
2) any bonding of AX ions with oxide ions would neutralize 
negative oxide ion charge as well as that of positive AX 
ions; 
3) with the above occurring, the process would slow down 
as the negative oxide ion layer is depleted (i.e. the 
intrinsic electric field diminishes) ; 


4) the negative oxide layer's influence on the surface 


oe 


unkwis to etotsond? SOws 
aongzstiib eivie eomie (VE asia See 
sstenists ae 9 etiam (8) 
ci jeiey od sneoe ShGSS le ymnMle Se 
gon sabioxtneio wa ads dealin Ae! aot Foatrt ts 
sit ida wmiididemes Ime maAee-biwow sa aaa 
iank start Sete Sit hanigeay, wiitankrq i ‘3a sit 
InaATU i ad fined “palbieser aren setort 6 
former twele Gok adt O39 ent “eLideanm: yis 
ibiue to 2ado0c0 ate ay hestsexd nail 
6.1) shine sr Yo eoaithetday een’ wid bluow yr 
'A-34 28 wt ane 46.4 eos ieeneta ah this 
'6 Slat? szrseels. an 4gn9e7. anti Bree teeter 
Howe to eer yet oct OVD Bly sande) pitied 
rioers aif o: (ott pret ae So 
C.D no Bm ts? | 
ies e679) et sees stil ssi lhctea 9 nee a 
| Repbe! fereveR “nt 
q lo Jneneror ee 19% egret « ebivorg’ a 
be weeirvved edd HE: 


vee 
ce 


Sct Lexwyen Bice cnot sites iw e1ot A to pntbaodt ae x4 
LA avittecg to tscy ee Liew Ge eerede nok in nt 


aweh #ole fives edacesy aes itniaind sil ee vy 
ony .9..i) Sotaiqoh staat nol ebixo avisspert ads es : 

| otfei sims bisks nivibelw\itantxint 

sistwvd ott no sehewkend e'aeval shixo avitapen orig’ (9 


3 | 


' : } ey ta ilar” a 

ig ca ten ee 
7 e » fn ] — a 7 ' 

ia we : be 7 


Hea: ‘ae 2 aaeliy We - 4 Dike 


f i A 


rs 


BiG Oral 
A model for the aluminum oxide barrier with the inclusion 


of a negative oxide ion layer is shown. 
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hydroxyls may serve to attract the positive hydrogen end 
of the hydroxyl group thus either canting it and/or 
hydrogen bonding with it somewhat and thus hindering its 
freedom of motion; 
5) any neutralization of the negative oxide layer by 
diffused AX ions will reduce the negative charge influence 
on the hydroxyl groups, thus allowing for reorientational 
motion away from the oxide layer toward the lead electrode. 
If an external field of forward polarity (AL posi- 
tive, Pb negative) is applied: 
1) the creation and migration of A& ions into the barrier 
is enhanced; 
2) the influence of the negative oxide layer on the surface 
hydroxyl groups is lessened due to a negative force 
attracting the hydrogen from the lead electrode side. 
Onwthe-other, handysifsaniieldvel oppesitespolarity:is 
applied: 
1) a reverse polarity field would compete with the nega- 
tive oxide ion layer for positive AX ions. Since a 
relative equilibrium existed prior to RP field application, 
the RP field would draw some already existing A ions out 
of the barrier; 
2) the positive force at the lead electrode, would hinder 


reorientation of surface hydroxyl groups. 


6.6. Annealing of AL-A% Oxide-Metal Junctions 


Attention is now turned to the influence of the top 
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metal electrode on the annealing process. In this experi- 
ment one junction of a pair of AL-AX oxide-metal junctions 
has a lead cover layer and serves as a control, while the 
other junction has a top electrode of a different metal, 
suebeas Ag, (AU, Au, Cu, In, Mg or Sn. ‘These junctions are 
thermally annealed. As will be seen, barrier parameters and 
annealing behaviour are often quite different from those 
observed when a lead top electrode is used. 

Fig. 6.22 shows a plot of the resistance ratio 
Ri/Ro VEUESUS ERG Octal, annealing time <t) for junctions wren 
various cover electrodes where Ry is the resistance at 
time t. The resistance behaviour of Pb, In, Sn, and Mg 
appear to be fairly similar. Cu may be included also in 
this group although its resistance fluctuated up and down 
with no set pattern - it also tended to develop filaments, 
with a consequent resistance drop, which were "burnt out" 
when the sweep bias was applied to the junction at helium 
temperature resulting in the resistance increasing to the 
"expected" value. The A& and Ag behaviour is much more 
dramatic than that of Pb while Au is in between. The dis- 
continuity sat) ~ 95 hours in the’ geldsjunction may be due 
to the “burning out” of a, filament. sai tnougheenesAY 
resistance ratio increase is much larger than that of Pb, 
the curvature of the graph is the same (i.e. an increase 
which tends to eventually saturate). The silver electrode 
VuNnctLoOn Ss curvature as radically different - it is reversed 


and shows no sign of resistance saturation even after 300 
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Pige orice 
Plots of RL/R, vs t for several AX-A& oxide-metal 


junctions, where R, is the resistance at time t, are 


uid 
shown. Data for Mg and Cu aresnot plotted (for the 
sake of clarity in the figure) but lie in the region 


bounded by the Pb and Sn curves. 
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hours of annealing. 

Since these are preliminary results with data for 
only one or two specimens of each cover electrode avail- 
able, interpretation of results is difficult at this 
stage, however, several observations can be made. 

During the preparation of these samples, the top 
electrodes were all initially evaporated very slowly so as 
tow gently cover the insulating barrier. In view of the 
Pease work (51, 92) on ionic radii of metals and depth 
Penetration into the barrier, a plot of normalized (to the 
lead control) 450 meV peak intensity versus ionic radius 
for metals in selected oxidation states is shown in Fig. 
6.23. The 450 meV peak is chosen because it is due to O-H 
stretching of hydroxyl groups on the surface of the 
aluminum oxide. Since these hydroxyl groups are proximate 
to the cover electrode, penetration of the cover electrode 
metal into this hydroxyl layer is expected to shunt out 
some of these groups and hence reduce the intensity 
observed in the IET spectrum. 

Previous work (92) on the use of ethylene (C5Hy) 
plasma discharge barriers used the 360 meV peak due to C-H 
stretching, while an even earlier study (51) on oxide 
barriers again used the 360 meV peak, this time arising 
from hydrocarbon contamination on the oxide surface. 


These works showed a correlation between the normalized 


360 meV peak intensity and ionic radius although the chosen 
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A graph of normalized 450 meV peak intensity 
Fy (450) /F pp (450) versus r, the ionic radius of the metal, 
is shown. Here Fy, (450) is the 450 meV peak intensity 


for a junction (with) cover metal. 
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oxidation states of the metals were not the same. This 
cursory study using the 450 meV peak also shows a 
relationship between the peak intensity and ionic radius. 
The ionic radii used in this and previous works are shown 
Tae taple. 6:1: 

Another way of looking at the penetration is by 
comparing the resistance of a junction with a cover 
Payer other than lead to the lead control. If there is 
more or less penetration, the resistance should be 
lower or higher as the case may be. This ratio is 
Plotted) insFig. 6.24 using ithe same ionic radii’ as disitted 
Bpove for Fig, G@.23. ° There is a good correlation for 
the five metals Mg, In, Cu, Pb, and Au. A& and Sn are 
above the line while Ag is below. 

Consider first the case of aluminum. Its IET 
spectrum shows a hint of the 450 meV peak (Fig. 6.25), 
although, having a smaller ionic radius than Mg (which in 
Fig. 6.26 shows no 450 meV peak at all) it is expected 
to show none. Aluminum thus does not appear to penetrate 
as far as expected which may be due to the following 
reason. Some of the aluminum that does penetrate may 
anteract. with the hydroxyl groups or even the aluminum 
Yoxide, layer itself. Since aluminum 1s fairly weactive 
(which is why it is easy to make tunnel junctions with it 
AS) aebase layer), this possibility seems Itkely when 


considering the large resistance increases exhibited by 


A£-AX oxide-A& junctions as well as a possible availability 
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Fig. 6.24 
This figure shows a relationship between Ryan’ pp pp 
and r, the Lonic) radius of) the mecar. Ru is the 
initial resistance of a junction with top electrode 


metal M and Ay its area. 
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Fig. 6.25 


The IET spectrum of an AL-A& oxide-A& junction is shown. 
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The IET spectrum of an AL-A& oxide-Mg junction is shown. 


Note the complete absence of the 450 meV peak. 
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OL negative oxide ions thought to be responsible for the 
PumMetIon! Ss intrinsic electric field. “During the ‘evapora-— 
tion of the aluminum cover electrode, the small ionic 
radius aluminum may penetrate so far as to come in 

contact with some of the oxide ion layer. A reaction with 
a portion of this layer would create more aluminum oxide 
insulator and thus increase the junction resistance to a 
value larger than expected if the aluminum remained as 
"inert" metal ions or particles. The apparent penetration 
is thus not as great as expected. Also, it is likely that 
aluminum ions from the top electrode (once the junction is 
completed) are closer to and therefore more strongly 
attracted to the negative oxide ion layer (responsible for 
the junction's intrinsic electric field) than those from 
the base electrode. Hence, resistance increases during 
annealing due to AX ions reacting with and neutralizing 
the oxide ion layer, normally due to base electrode A& ions 
in AL-AX oxide-Pb junctions may well, in this case, be due 
to top electrode AX ions. As this happens, the oxide ion 
electric field is decreased and AX ions originally present 
at the base electrode will experience a less negative (i.e. 
more positive) potential forcing some of them out of the 
barrier. This will cause by to increase which is what 
happens (Fig. 6.27). The resistance decrease due to the 
loss of these ions from the barrier at the base electrode 
is small (for example, the application of a 075 V RP field 


to an A-AX oxide-Pb junction yields a resistance ratio 
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Fig, 1onec 7 
The behaviour of ry and 5 versus t for an AL-AX oxide-A 
Junction is shown. oy is represented by circles and 5 


by squares. 
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minimum <0 .95swith a change in aT comparable to that 
in the A-A& oxide-A& junction) compared to the overall 
increase due to the reaction occurring at the top electrode, 
therefore the net resistance change is one that increases. 

In the case of the tin cover layer, only one sample 
is available at present so no suggestion as to its behav- 
Tour will be, put fonward. 

Unlikevaluminum or tin,, gunctions with a silver top 
eeectrode are, below the dine (Fig. 6.24), that is to say, 
the silver penetrates further than expected. Its resistance 
behaviour (Fig. 6.22), as well as changes in 1 and 5 
(Fig. 6.28) are quite different from the other junctions. 
The rapid initial decrease in 5 implies that the silver 
Side of the junction is becoming more positive. This, 
coupled with the eventual resistance increase suggests 
that positive Ag ions are being created at the silver elec- 
trode and migrate into the aluminum oxide (perhaps under 
the influence of the intrinsic electric field). This 
influx of positive ions could cause the A electrode side 
to experience a positive potential thus forcing A ions 
Out of the barrier and causing rT 40 rise. ‘One notes 
however, that 5 does not saturate but keeps on climbing; 
gilso, the resistance Continues to Gise] 9G 1S Seen that 
5 reaches a minimum which is maintained for 100-150 hours 
and then increases again. This suggests that some of the 
Ag ions have moved into the oxide and eventually move far 


enough into the oxide so as to decrease positive ion 
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Fig. somec 
This graph shows how ory and 5 (barrier heights at the AX 
and Ag electrodes) change with time for an A&-A oxide-Ag 
junction. The Ag electrode is seen to become rapidly 


positive in the initial stages of thermal annealing. 
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influence on the oxide-Ag interface (thus causing it to 
become less positive). Meanwhile, as the silver ions move 
Rowardpthe Avselectrede ,<Av ions careicontinua My being 
forced out thereby causing the increase in i: 

The question one asks is why does this happen to 
Silver and not to the other metals? Perhaps the answer 
lies in the work on B-alumina (81). Beta-alumina is an 
aluminum oxide doped with a small amount of NaO (empirical 
foxumua Na 40 «LIAL <Os)) rGL2.:, cb8)ssandesomnamed ‘because, it was 


2 eae 


thought to be an isomorph of ALO Since the presence of 


5) 
sodium was either unknown or considered unimportant by 
early workers. “ The Na ions in this material, of interest 
rm theritveld soft fast wwonic ityansporviintsolids;, iare highly 
mobile. It is also relatively easy to substitute Ag ions 
for the Na ions in this material. Other ions used for 
substitution (which substitute less easily than Ag) tend 
to be monovalent cations (in particular, the alkali metals) 
unlike most of the metals used for the junction cover 
electrodes. In addition, work on 8-alumina shows that 
substitution of Ag for Na causes the lattice to contract; 
ali other cationic substitutions resule in a lettice 
expansion. With these points in mind it is interesting to 
speculate that while the aluminum oxide of the barrier is 
not likely to contain sodium and is thus not f-alumina 

(but rather an amorphous alumina (125)), it may possess 
properties necessary to induce Ag ions to enter into the 


oxide. A tendency of Ag to enter into the insulator could 
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account for its initial lowered resistance compared to Pb 
md 6.22)). During the first 100 hours of annealing a 
buildup of positive Ag ions may occur within the aluminum 
Oxide near the oxide-Ag interface causing a decrease in 5 
along with increases in resistance and barrier thickness 
ee(Pag. 6.29). Inthe next 100 hours, Ag ions may predomi- 
nantly drift into the aluminum oxide, causing the greater 
resistance increase, and enter into appropriate sites to 
effect a lattice contraction and hence a decrease in s. 
During this stage, 5 is faLciy constant indicating that 
any positive charge that drifts into the oxide is replaced 
by new ions. Following these first 200 hours, 5 starts 
to increase again, suggesting the number of Ag ions near 
the top electrode is decreasing as the ions enter into 

the oxide, meanwhile still increasing the resistance and 
decreasing s. Although one can thus interpret these 
annealing results, questions remain and further work needs 
to be done on A&-AL oxide-Ag junctions before making any 


definite conclusions. 
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The change in barrier thickness with time, s vs t, for 
AX-AX oxide-Ag is shown. At first, Ss increases to a 


maximum and then decreases quite rapidly. 
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CHAPTER 7 


SUMMARY AND CONCLUSIONS 


Veo CONCLUSIONS 

In the course of studying annealing effects in 
tunnel junctions, many interesting observations were 
made. 

AL-AX oxide-Pb tunnel junctions subjected to 
thermal annealing at room temperature show overall 
increases in resistance and in IET peak intensities (118 
and 450 meV peaks) with respect to time. A reorientation 
of surface hydroxyl groups under the influence of an 
intrinsic electric field of forward polarity has been 
suggested. 

Voltage annealing of AX-AX& oxide-Pb tunnel junctions 
tended to enhance or inhibit the changes observed in the 
thermal annealing case, depending upon the polarity of 
the applied voltage. Enhancement was observed with a 
forwards polarity.bias... In facies the application, of. a 
forward polarity voltage (0.5 V FP) causes an accumulation 
of positive charge, presumably Af ions, at the A&-AX oxide 
interface. On the other hand a reverse polarity bias 
(0.5 V RP) causes the removal of some positive charge 
from the AL-AX oxide interface region. In addition, peak 
intensities do not change because the magnitude of the 
bias seems sufficient to negate the junction's intrinsic 


field, thereby leaving no electric field to reorient 
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molecular groups. Junction aging is thus postponed until 
the reverse polarity voltage is removed, at which time 
normal thermal annealing occurs. 

Voltage annealing with alternating polarity showed 
that a large portion of the positive charge that could be 
moved into the aluminum oxide under the influence of a 
forward polarity bias remained trapped within the barrier, 
while a smaller portion was fairly mobile and could be 
moved in or out depending upon the polarity of the applied 
bias. The largest changes occurred within the first two 
hours of annealing. 

THe “Origin” or the™ junction’ s, rnhtrinsic erectric 
field, estimated to be 0.5 V FP, has been suggested as 
being due to a layer of negative oxide ions near the 
aluminum oxide surface, just below the surface hydroxyl 
groups. These oxide ions are assumed to be formed as part 
of the process of aluminum oxide growth in a plasma 
discharge. 

Preliminary thermal annealing studies of AL-AxX 
oxide-metal tunnel junctions resulted in some behaviour 
Simivar to Po tor Cu, in, Mg, ana on, wide AG, 4, ana 
Au were different. Speculative models were presented for 
the aluminum and silver behaviour, but further work needs 
to be done on these, and the other metal cover electrodes. 
Finally, a relation between the surface hydroxyl stretch 
peak intensity (450 meV) and ionic radius of the top 


electrode metal showed that more penetration into the oxide 
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occurred the smaller the ionic radius became. 


7.2 Afterthoughts 

At the outset, the A&-A2 oxide-Pb system was chosen 
as the system to study since it is the most common type 
of junction used in IETS, the oxide being frequently 
doped with organic molecules. The oxide barrier, being 
less complex than an oxide plus organic molecule barrier, 
was thought to be a "Simple" material to study. However, 
through the course of experiments described in this thesis, 
as well as research into the literature concerning areas 
such as catalysis, chemistry, mineralogy, IR and Raman 
spectroscopy, etc., the "Simple" barrier was found to be 
somewhat more complicated. 

Aluminum oxide, hydrated and unhydrated, exists in 
a multitude of forms, both man-made and natural (34, 35, 
49, 124, 125). Only one form of unhydrated aluminum oxide 
existsrindnature.™.Thistis ‘the most stable trorm, td=alumina 


(a-AX 03) which is also called corundum and is one of the 


2 
hardest known minerals. Hydrated forms of alumina (124) 
found in nature include boehmite and diaspore (AXO(OH)), 
and gibbsite and bayerite (AR (OH) 5). Heating and 
dehydrating these minerals under various conditions 
yields aivariety of other aluminas’ including y7 "67" n; 
pe yor, Vandstinally;¢withtenoughsheat a-alumina is 


arrived at (35). y-alumina is of interest as a catalytic 


substrate (88, 89) and has been well studied (103, 104). 
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The alumina formed by the electrolytic anodization of 
aluminum metal has been found to be amorphous alumina 
(125) which can undergo a transformation to y-alumina 
upon heating above ~ 400°C. The aluminum oxide formed in 
the plasma discharge is most likely amorphous alumina 
(125) covered with some surface hydroxyl groups. 

Infrared spectroscopy has been used to study the 
Various fforms) oftaluminia 1038, 45 88=9.0:cel03=105) 120-123, 
Eos) MA iBandsiloccur Hinethe 30-120) meV region ofthe various 
spectra for different types of alumina. Other work has 
also shown that the absorption region for metal-oxygen 
stretching (in our case A£-O stretching) depends upon the 
coordination of the metal with respect to the oxygen (ll, 
45,0 203) Av For itetrahedral coordination (ALO, ) the region 
is 85-110 meV (700-900 em7+) whereas octahedral coordina- 


tion (A20,) has the range 60-80 meV (500-650 em™+).. 


2) 

An X-ray diffraction study on amorphous alumina 
(a-alumina) by A.M. Jones (Ph. D. thesis, 1974, Texas 
Christian University) suggests that most of the aluminum 
is in tetrahedral coordination. Also, as mentioned in 
Chapter 5.2, work on anodized alumina (an amorphous 
material) gives infrared bands at 960 and 650 om > (127 
and 79 meV) due to alumina modes. 

The identification of the 118 meV peak as being 
due to A£-O stretching is thus reaffirmed, however, one 


now again must consider the 85 meV shoulder as being 


composed of modes other than OH bending, such as modes 
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due to tetrahedrally coordinated and some octahedrally 
coordinated alumina. The resolution of this will 
certainly require more study and perhaps an examination 
of the infrared spectroscopic behaviour of other metals 


Eenvcetrahedral@coordinationn (27.,,0290! 99; mMigs 121.2128) 


7.3. Suggestions for Further Work 

In' the course of this’ work a better insight into 
the nature of the aluminum oxide barrier in tunnel 
junctions has been gained. However, as is usual, more 
questions have been raised than answered and so a few 
suggestions for further work will be presented. 

As mentioned in Chapter 6.4 on alternating polarity 
voltage annealing further study should examine 1) the 
length of the annealing period ay 2) the magnitude of 
the annealing voltage, and 3) the effect of initially 
annealing the junction with a reverse polarity instead of 
forward polarity bias. Of interest would also be the 
application of "asymmetric" voltages, such that the net 
field (sum of the applied andwinmtrinsic fields) would be 
the same in FP and RP annealing. 

On the whole, more work needs to be done on 
junctions with different cover electrodes. Junctions with 
Ag and Ag top electrodes are of particular interest in 
that their annealing behaviour is markedly different from 
Pb. As a note, the IET spectra of Ag junctions suggest 


the existence of an undershoot on the 450 meV peak 
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(surface OH stretching). /This’may be related to the 
predicted undershoot of Davis (33), and Birkner and 
schattke=(I6) ; 

It has been suggested here that the 85 meV 
shoulder is due to more than OH bending modes. Infrared 
work suggests the presence of other alumina modes in 
this region. A study of the annealing behaviour of 
AX-AX oxide-Mg and AX-AX oxide-AX junctions may provide 
the first step in separating OH bending modes from 
alumina modes, since the small ionic radii of these 
metals causes the surface OH groups to be shunted out. 

Finally, one eventually hopes to study metal- 
insulator-metal junctions with insulators other than 
aluminum oxide in order to better understand the inter- 
action of the different metal electrodes with the 


insulating barrier. 
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| M. K. Konkin and J. G. Adler 


| PACS numbers: 73.40.Gk, 74.50.+r, 85.25.+k 


Recent interest in inelastic electron tunneling spec- 
troscopy (IETS) has necessitated quick and convenient 
fabrication techniques for stable tunnel junctions. How- 
‘ever, in the study of junctions from 0 to 500 meV, in- 
Stabilities are often seen. The instability shows up as a 
large noisy patch in the second derivative signal d’J/dV’ 
=do/dV, where / is the current, V is the voltage, and 
0 is the conductance. Figure 1 shows three tunneling 
Spectra with instabilities present. 


The junctions in this study were prepared by first 
evaporating an aluminum base film onto a clean glass 
Slide. A plasma discharge of either oxygen, methane, 
or oxygen followed by methane formed the barrier. 
‘Although this work deals only with three types of plas- 
’ma-prepared barriers, these instabilities are also seen 
‘in junctions with thermally grown, liquid-doped, and 
vapor -doped barriers. After evaporating lead top 
electrodes, the junctions were cooled to 77 K in liquid 
nitrogen. The spectra for these junctions were obtained 
lat 4.2 K over a range 0—500 meV using an updated 
version of the bridge-minicomputer system described 
by Adler and Straus. ! 


The data collected often showed instabilities. By 


*Supported in part by the National Research Council of Canada, 
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Instabilities in thin tunnel junctions?) 


Department of Physics, University of Alberta, Edmonton, Alberta T6G 2J1 Canada 
(Received 21 November 1977; accepted for publication 22 January 1978) 


Tunnel junctions prepared for inelastic electron tunneling spectroscopy are often plagued by instabilities in 
the 0-500-meV range. This paper relates the bias at which the instability occurs to the barrier thickness. 


looking at several dozen junctions it was found that the 
position of the instability was related to the product RA 
(which is proportional to the barrier thickness’) where 


E (cm*!) 


ce) 1000 2000 3000 4000 


(meV!) 


| do 
o> dV 


103 


0) 100 200 300 400 500 


E (meV) 


FIG, 1, These three offset curves (Al negative) illustrate the 
shift of instability position with increasing RA, For these junc- 
tions which were prepared by putting different amounts of 
methane on an oxide, RA is (a) 1.83 2mm’, (b) 4.46 Qmm?, 
and (ce) 8,23 2mm, 
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IG, 2, This figure shows the instability position, V; (center 

f the instability), plotted on a quadratic scale against RA, The 
jircles represent junctions with an oxide barrier, while the 
quares are for barriers made of methane plasma discharged 
Wa an oxide, Note that the methane plus oxide plot will inter- 
ept the ordinate at a positive value, while the single layer 
|»xide) appears to pass through the origin, 


jalues of RA were obtained by varying the oxygen and 
jrethane plasma discharge times during junction fabri- 


y position and RA suggests that the explanation of 
ese instabilities involves the power dissipation of the 


_ All the data shown were from junctions studied within 
|day or two of preparation. There are indications that 
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the instability position changes with thermal cycling 
and aging, shifting it to higher bias. 


When an instability is present it seems to affect the 
IET spectrum only as far as noise is concerned. At 
biases below the instability the signal is relatively 
quiet, while above it the noise level is increased 
(Fig. 1). The widths of observed instabilities range 
from a few meV up to ~100 meV. Another property of 
the instability is that it moves to lower and lower bias 
voltage as the magnetic field applied to a junction is 
increased, vanishing near or above the critical field of 
bulk Pb (e.g., ~1 kOe for an instability of ~ 280 meV 
zero field position). 


In conclusion, the thinner the barrier the lower the 
bias at which instabilities occur, Because this often 
masks peaks of interest to IETS, junction preparation 
parameters should be adjusted to produce a thickness 
(i.e. , an RA value) such that instabilities will not be 
present in the region of interest. If RA is chosen so that 
instabilities occur only above the region of interest, 
reasonable signal-to-noise ratios are obtained, so long 
as junction resistances are kept low. Finally, the 
quadratic relationship between the instability position 
and RA suggests that an instability model should involve 
power dissipation in the junction. 


1J.G, Adler and J, Straus, Rev. Sci. Instrum, 46, 158 
(1975). 

’R. Magno, M.K, Konkin, and J.G, Adler, Surf. Sci. 69, 
437 (1977). 

‘It should be noted that the center of the instability V, differs 
slightly depending on bias polarity with V, for Al positive 
being closer to zero bias by about 10 meV or less. 
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A STUDY OF INSTABILITIES IN SUPERCONDUCTING TUNNEL JUNCTIONS*+ 


M.K. Konkin and J.G. Adler 


Department of Phystes, Untverstty of Alberta Edmonton, Alberta, Canada T6G 2J1 


Résumé.- Les jonctions 4 effet tunnel supraconducteurs présentent souvent des 
instabilités. Cet article relie la tension 4 laquelle l'instabilité se produit 
a 1'épaisseur de la barriére et démontre la dépendance de la position de l'ins- 


tabilité au champ magnétique. 


Abstract.- Superconducting tunnel junctions often exhibit instabilities. This 
paper relates the bias at which the instability occurs to the barrier thickness, 
as well as showing the magnetic field dependence of the instability position. 


Superconducting tunnel junctions find 
many applications in the study of : phonon spectra 
of superconductors, inelastic excitations in the 
junction barrier, etc...It is therefore of prime 
importance to have quick and convenient fabrica- 
tion techniques for stable tunnel junctions. An 
instability shows up as a large noisy patch in 
the second derivative = = = where I is the 
current, V the voltage, and o the junction conduc- 
tance. The insert in Figure | shows three tunnel- 


ling spectra with instabilities present. 
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Fig. 1 : This figure shows the instability position 
Vi (center of the instability), plotted on a qua- 
dratic scale against RA. Circles represent junc- 
tions with an oxide barrier, while squares are for 
barriers made of methane plasma discharged on an 
oxide. The three offset curves (Af-ve) in the 
insert illustrate the shift of instability posi- 
tion to higher bias with increasing RA. 


——— el 


+Supported in part by the National Research Council 
of Canada. 


The junctions in this study were prepared 
by first evaporating an aluminium base film onto 
a clean glass slide. A plasma discharge of either 
oxygen, methane, or oxigen followed by methane 
formed the barrier. Although this work deals only 
with three types of plasma prepared barriers, 
these instabilities are also seen in junctions 
with thermally grown, liquid-doped and vapor- 
doped barriers. After evaporating lead top elec- 
trodes, the junctions were cooled to 77 K in li- 
quid nitrogen. The spectra for these junctions 
were obtained at 4.2 K over a range of 0-500 meV 
using an updated version of the bridge-minicompu- 
ter system described by Adler and Straus /1/. 

The data collected often showed instabili- 
ties. By looking at several dozen junctions it was 
found that the position of the instability was re- 
lated to the product RA (which is proportional to 
the barrier thickness /2/)where R is the zero bias 
resistance at 77 K and A is the junction area. 
Figure | illustrates how the instability position, 


V., shifts with a change in RA3. At a high enough 


’ 
ae of RA one no longer sees instabilities in 
the 0-500 meV range. Figure 1 also shows the ins- 
tability position, V;.0n a quadratic scale against 
RA. The different values of RA were obtained by 
varying the oxigen and methane plasma discharge 
times during junction fabrication. The quadratic 
relationship between the instability position and 
RA suggests that the explanation of these instabi- 
lities involves the power dissipation of the junc- 
tion. A study of the behaviour of the instability 
in a magnetic field (Figure 2) shows that it moves 
to a lower and lower voltage as the field is in- 


creased. The plot is shown in reduced field units 


39 


tus 
(PL 


ue Ae): TE ie 


satel Sasveeg?e 2. poe 
sigeng oe 09s idk 
aa ey warsreou oS: ay ; 


Py ad Parie tee 2 TP it 
Seas itt Fobvaae bt} 92 ypey tee 
rf’ ieoe Ve i) Yeadon “rs w i 


meoe orad Yau a hd L @peed | Src), alt. 

Oyu we \eute fe . GrDtweggers pared ge 
|, ring & Su he Walp e 
our Sow wit Bepadto? msgsen 7 @ nbdicer . ep ‘ 


. a 20) 
a a: me im 
snow ates cepa Sh ce. Tee we 


vhrviad bayegesy apehlg >) KAY) Sears abl os 


a 
nh} ya L Ohne opts OTe sca apa devity, : Tied ; 
pon ko tke Leteeedod el are hanes : rm 
» yam tes fr eae asF02 ee bi 7: of 4 a : say a 
ch ak MN AF BETSRS aaa eos hi wes | 
ufietsceep been “MOA nig se sat om muy Je ae ; 
Vals 007+) Selapake a abe F32F 1 Sem : iW be a 7 es 
“yeahs i+ ie-sRbhed’ 48420 Ae aes ba telgy a eae Teva oe 
is i ‘erag dd en PeeOe. ee Ge, ee ee ' | = ie 
-glédusen! bewaly peeto tasass foo By 1G, +t : 


‘ hr 
2au'3h CUGlIswl haa feieree ae a i vA. sit , are Raa aye 
ise ‘w é 

=) a 


: 
id 


= 


ey Raw ¥ IF) ides pis wit) in bos theog # 


wd Legs saad que B27 apdoet) AG erqeng ea. ot Pra 
aid OCs ate 848 p parts 484 parasol a? on 
een oie ott 23 bbe a aN 1A aap mm. 
ohS be uty v2 bad aul gad Nae ee: veel ie 
dpucae dnote a8 ty G7 agealls de a adzelien gy 
of épful liseians' sees) Seg7ol on i ae 
~¢pd ets, Gta jath { avvagtt . aba ae 
sonjeat ahem 2e4nubhp oe 80, 4% tees este 
yt bettsr aun: eee aia tu: eer: sede 
epee sh mapas wrgdige imal oO : 
dis@yeaup bet asi eat tots sn girl 
fee nerd cee? coi ganar ot, had 
“laren esa? to galsaceiahe ot yess 
-seyt =i) te peagaieteniel mits ny 


209 
C6-592 


to take into account the slight variation in criti- References 
cal field from junction to junction. 


/\/ Adler, J.G. and Straus, J., Rev. Sci. Instr. 
46 (1975) 158. 


/2/ Magno, R., Konkin, M.K., and Adler, J.G., Sur- 
facesSsei. 69 (1977) 437. 


/3/ It should be noted that the center of the ins- 
tability V; differs slightly depending on bias 
polarity with Vz for A& positive being closer 
to zero bias by about |]OmeV or less. 


V4 (H)/V_(0) 


0 2 A 6 8 1.0 
H/He 


Fig. 2: Instability position vs applied magnetic 
field for two different oxide junctions is shown 
here in reduced units. V{(H) is the field dependent 
instability position, Vyz(0) is the zero field va- 
lue - for circles VI(O) = 173 meV while for squa- 
res it is 279 meV - and A. is the critical field. 


Instabilities present in a junction seem 
to affect the spectrum only insofar as noise is 
concerned. In general the noise level is low at 
voltages below the instability while being some- 
what larger above it (see insert Figure 1). The 
widths of the instabilities have been observed to 
range from about ImeV to about lOmeV. The data 
presented here was obtained within a day or two of 
junction prepatation. There is some evidence that 
the instability position shifts to higher bias 
with thermal cycling and junction aging. 

In conclusion we have shown that there is 
a quadratic relationship between the instability 
position and RA, indicating that power dissipation 
in the junction is involved. The magnetic field 
dependence of the instability position, as well as 
the absence of instabilities in junctions with 
both electrodes normal indicates that superconduc- 


tivity is involved. 
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EFFECT OF COVER ELECTRODE METAL ON INELASTIC ELECTRON 
TUNNELING STRUCTURE * 


R. MAGNO, M.K. KONKIN and J.G. ADLER 
Department of Physics, University of Alberta, Edmonton, Alberta, Canada T6G 2J] 


Received 13 June 1977; manuscript received in final form 8 August 1977 


Inelastic electron tunneling spectra peak intensities have been measured for a group of 
AI—M,, junctions for which My was Ag, Au, In, Pb, Sn or Zn and the barriers were formed in an 
ethylene glow discharge. The intensity of the 360 meV peak which is associated with the 
stretching of carbon—hydrogen bonds was found to increase with junction barrier thickness for 
a given cover electrode. The intensity decreases as the cover electrode is changed from Ag to Pb, 
Sn, Au, In, and finally Zn. This order correlates with decreasing ionic radius, suggesting that the 
decrease in peak intensity is due to the increased penetration of the cover electrode into the 
insulator, resulting in a thinner barrier. 


1. Introduction 


It has long been known [1—S] that Pb is the optimum metal for use as a cover 
electrode in metal—insulator—metal tunnel junctions prepared for the purpose of 
studying inelastic electron tunneling. This paper deals with the variation in intensity 
of inelastic electron tunneling spectra with junction counter electrode. Data are 
shown for a series of Al—M, junctions where M, is either Ag, Au, In, Pb, Sn or Zn, 
and the barriers are made in an ethylene glow discharge [6]. Geiger, Chandrasekhar, 
and Adler [7] examined the intensity of inelastic structure for various cover elec- 
trodes using oxygen glow discharged aluminum base layers and showed a correla- 
tion between the intensity and the ionic radius of the cover electrode. Earlier work 
of Handy [8] which studied resistances as a function of the covering electrode for 
AI—M, junctions indicated a correlation between the atomic radius of the counter 
electrode and the resistance of the junction. 

Inelastic electron tunneling manifests itself by an increase in conductance, o, at 


energies |ev| 2 hw, where fiw is the threshold energy for the inelastic effect in ques- 


tion [1—5]. The conductance increases since each inelastic excitation adds a new 
channel to the tunneling process. In the case of a band of molecular vibrations the 
change in conductance, 60, is rarely more than 2%. For this reason it is advantage- 


* Work supported in part by the National Research Council of Canada. 
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ous to measure the derivative of the conductance with respect to energy do/dv. In 


this work we have used a measurement system similar to that described by Adler 
and Straus [9] to obtain values of the intensity 


1 do 
=f [2 = at] du, (1) 


where 09 is a value of the conductance at low bias chosen for normalization and 
g(v) is the smooth background do/dv under the inelastic peak. 


E (meV) 


OQ Al-Ethylene-Pb 
J1846 


5 500 
E (meV) 


Fig. 1. Conductance (a) and its derivative (b) for a typical Al—Pb junction. og is the 
value of o at 35 meV and is used to normalize the curves. 
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2. Experimental 


The samples examined in this study consisted of a pair of junctions formed on a 
common base layer. The base was made by evaporating an aluminum film onto a 
clean glass slide. The barrier was produced in a dc glow discharge of ethylene. 
Details of the glow discharge procedure have been reported elsewhere [6]. The two 
junctions sharing a common base layer and insulating barrier were completed by 
evaporating a Pb control electrode for the first junction and a metal M, for the 
second junction. Measurements were carried out at 4.2 K in liquid helium and the 
data were taken on an improved version of the bridge and minicomputer system 
described earlier [9]. 


14 Al-Ethylene-Pb 
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Fig. 2. (a) (1/og) do/dv near 360 meV witn the Pb electrode biased positively for an 
Al—Pb junction. The dashed curve is the background g(v). (b) (1/09) do/dv with 


background removed, this curve forms the integrand of eq. (1). 
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The junctions were sewpt from 0—5S00 meV (M, positive) except for those with 
a Pb cover layer for which the region below 35 meV was omitted to avoid satura- 
tion of the lock-in amplifiers by structure due to the Pb phonons and the supercon- 
ducting energy gap. 

A typical conductance curve is shown in fig. 1a, while its derivative (1/09) da/du 
is illustrated in fig. 1b. These data and all those shown in this paper have been 
normalized at 35 meV to allow for comparison of various junctions regardless of 
junction resistance (the conductance at 35 meV differs from that of zero bias by 
less than 3% in all cases). A typical high resolution sweep in the 275 to 425 meV 
range (Pb electrode positive) is illustrated in fig. 2a, along with the dashed back- 
ground curve, g(v). The integrand as used in eq. (1) is shown in fig. 2b. The mini- 
computer based system [9] allows this background subtraction and the ensuing 
integration for obtaining intensities to be carried out in a few minutes. All junction 
pairs were analysed in this way to determine the intensity of the 360 meV peak. 


3. Results and discussion 


Consider first the Al—Pb junctions. The resistance of a junction should increase 
exponentially with the thickness of the barrier, thus a plot of the intensity versus 
In(10RA) (where R is the junction’s resistance and A its cross-sectional area) may 
be expected to be linear. Such a plot is illustrated in fig. 3. It should be noted that 
the junction areas varied by almost an order of magnitude from less than 0.1 mm? 
to greater than 0.7 mm?. The results thus indicate that the tunneling measurement 
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Fig. 3. Intensity of the 360 meV peak plotted against In(10RA) for a group of 
Al—Pb junctions. 
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Fig. 4. Typical data for the various covering metals used. 


is sampling carbon—hydrogen bonds throughout the volume of the barrier rather 


than just those near a metal—barrier interface. 
We turn our attention now to the variation of intensity due to the cover elec- 


trode metal. Those results are shown in fig. 4 and plotted against ionic radius in fig. 5. 
These results suggest that smaller ionic radius metals penetrate the barrier more, 
effectively resulting in a thinner barrier junction. 
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Fig. 5. Intensity of the 360 meV peak versus the ionic radius of the cover metal. 
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4. Summary 


In conclusion, we found that the intensity of the spectra increases with barrier 
thickness for a given cover electrode (fig. 3), while it depends on the penetration of 
the cover electrode into the barrier when different covers are used (fig. 5). The 
distance the metal penetrates depends upon its ionic radius. This conclusion is sim- 
ilar to that reached by Geiger et al., except that the data (fig. 5) are plotted using 
ionic radii different from those previously used. A comparison is given in table 1. It 
should be noted that the intensities measured here are those defined by eq. (1), 
while the data of Geiger, Chandrasekhar, and Adler [7] refer to the relative inten- 
sities of the hydrocarbon peak near 360 meV to the alumina hydrate (or alumina) 
peak near 118 meV. The ionic radii used in the earlier work [7] were those given by 
Pauling (which list only one oxidation state per element), while we have chosen 
more frequently occurring oxidation states with perhaps the exception of tin. If 
one replots the data of Geiger et al. one finds that the radii used in this work 
actually give better agreement. The important point of both papers is a qualitative 
one which indicates that the ionic radius is an important parameter in determining 
the properties of tunnel junctions. 

The details of the penetration of the cover metal into the barrier merit further 
examination. At this point it is not known whether the penetration occurs as the 
cover metals are deposited or by subsequent diffusion. 

It is possible that this experimental technique could be developed further as a 
method for determining the ionization state of metallic ions in the type of envi- 
ronment existing in tunnel junctions. 

Finally it should be pointed out that from fig. 5 one might expect silver to be 
the optimum cover metal; this however is not so since Al—Ag junctions are noisier 
than Al—Pb junctions as well as having steeper conductances. 


Table 1 
Summary of ionization states and ionic radii 
This work Geiger et al. [7] 
Radius (A) [10] Radius (A) 
Agt! 1.26 Agt! 1.26 
Pot? 1.20 Bone 0.84 
Sn? tune Sa 0.71 
Au*3 0.85 Au?! 1.37 
In’? 0.81 In 0.81 


Zn*2 0.74 
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Experiments show that for an organic material placed on an oxide one sees 
an asymmetry of inelastic electron tunneling peak intensity for bias volt- 
ages of opposite sign. When the barrier is composed only of the organic 
material with no deliberate oxidation such asymmetry disappears. 


RECENTLY there has been considerable interest 
regarding the asymmetry of IETS (inelastic electron 
tunneling spectroscopy) peak intensities [1]. This 
asymmetry is manifest in the unequal intensities (arising 
from molecular vibrations of organic molecules in Al— 
oxide—organic—Pb junctions) observed at opposite 
junction bias polarity. When the Al base film is negative 
with respect to the Pb top film one obtains a larger peak 
intensity than with reverse polarity. 

Yanson et al. [2] first looked at peak intensity 
asymmetries in Sn—Pb and Al—Pb junctions with an 
organically doped oxide barrier. They showed that an 
asymmetry does exist and gave reasons to support this. 
Since they did not use calibrated data (i.e. the second 
derivative do/dV, where a is the junction conductance 
as opposed to the 2nd harmonic voltage signal) little 
more on the subject could reliably be said. 

Kirtley, Scalapino and Hansma [3] did theoretical 
calculations for peak asymmetries. To support their 
theory, ratios of conductance change Ao(+ V)/Ao(— V) 
vs voltage V were found. The Ao(V) refers to the height 
of a step in the conductance curve caused by the opening 
up of a new tunneling channel due to the inelastic 
excitation of a molecular vibration, for example. For 
Al—oxide—benzoic acid—Pb data, a plot of Ao(+ V)/ 
Ao(— V) vs V was shown to lie near their theoretical 
curve. 

In this work we use absolute IETS peak intensities 
obtained from calibrated second derivative data to show 
the asymmetry. Moreover, we demonstate the lack of 
asymmetry for a barrier composed of only the organic 
material with no prior oxidation step. 

Tunnel junctions were prepared by first evaporating 
an aluminum base film onto a clean glass slide. The alu- 
minum film was then oxidized in an oxygen—helium 
plasma discharge followed by a methane discharge to 
cover the oxide with a simple organic material. Finally, 
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a lead top electrode was evaporated to complete the 
junction (Al—Pb was chosen for best resolution) [4]. 
The junctions thus produced were studied at 4.2 K and 
the data were obtained using an updated version of the 
bridge-minicomputer system described by Adler and 
Straus [5]. Since we calibrate our data we are able to 
work with absolute quantities such as do/dV = d7//dV? 
as opposed to raw second harmonic signals, or relative 
units normalized to some arbitrary peak height. 

For the asymmetry study we chose the second 
derivative peak due to C—H stretching modes occurring 
near 360 meV because of: (1) its prominence, (2) iso- 
lation from other peaks, and (3) the relative ease with 
which a background can be subtracted. 

The peak intensity is calculated from the data by: 


where do is the conductance used to normalize da/dV 
(in this paper 09 = o(— 30 meV), since we did not sweep 
through zero bias in order to avoid the superconducting 
energy gap and large Pb phonon peaks) and g(V) is the 
smooth background under the inelastic peak [6]. We 
define [4] an asymmetry parameter a(V) as: 


CK A CAG) 


where F(+ V) is the peak intensity when the base film is 
positive with respect to the top film (i.e. for these junc- 
tions Al positive, Pb negative) and F(— V) is the peak 
intensity with the base film negative with respect to the 
top film (i.e. Al negative, Pb positive). When F(— V) is 
larger than F(+ V) (as typically occurs for organic 
vibrations in Al—oxide—organic—Pb junctions) a(V) <1. 

For junctions with plasma discharged methane on 
aluminum oxide, the asymmetry a(360 meV) ranged 
from 0.6 to 0.8 (several dozen junctions with varying 
oxidation and methanation times were studied). Junc- 
tions with no prior oxidation discharge gave a(360 meV) 
= ] showing that very little or no asymmetry exists for 
a barrier made of a single simple material. We have 
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Fig. 1. (a) do,/dV for an Al—oxide —methane—Pb 
junction. (b) do,/dV for an Al—methane—Pb junction. 
0, is the even part of the junction conductance. 


observed similar asymmetries with several formic acid 
vapor doped on plasma oxidized aluminum junctions, 
while formic acid vapor deposited on pure aluminum 
shows a near unity [7]. Typical results are shown in 
Figs. 1 and 2. 

These results can be understood as follows. In the 
asymmetric junctions (which have prepared oxide layers) 
the organic material is closer to the Pb electrode than to 
the Al one. When electrons tunnel from Al to Pb they 
cross most of the barrier at a high energy, then lose 
energy to an inelastic excitation, and have to tunnel 
only a short distance at low energy. Going from Pb to 
Al the electron loses most of its energy near the begin- 
ning of its journey and thus has to tunnel through most 
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PAES) 
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E (meV) 


Fig. 2. Peak asymmetry illustrated by plots of (1/09) x 
(do/dV) — g(V) (for the two bias directions) used to 
calculate the intensities F(+ V) and F(— V) as described 
in the text for the two junctions of Fig. 1. For (a) 
a(360 meV) ~ 0.81 and for (b) a(360 meV) ~ 0.98. 


of the barrier at low energy. Because the probability of 
tunneling diminishes as the electron energy decreases, 
fewer electrons will get through the barrier in the latter 
case. Hence, one sees a greater inelastic tunneling inten- 
sity for Al negative with respect to Pb than vice versa. 

In conclusion we have shown that by using absolute 
inelastic electron tunneling peak intensities one can 
study peak asymmetry. For barriers composed of an 
organic material on an oxide one sees a definite 
asymmetry, while for a barrier made of the organic sub- 
stance alone the asymmetry disappears. Our observations 
thus clarify the earlier confusion [1] which motivated 
this study. Finally it should be noted that this method 
may be used to determine whether a barrier is homogen- 
eous (a ~ 1) or heterogeneous (a # 1). 
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This paper deals with barrier parameters calculated from the coefficients 
of polynomial fits to electron tunneling conductance curves for a large 
sample of junctions with various types of barriers. Methods of data 
handling, along with a comparison of the results from two different para- 
meter calculations are shown. Finally, relationships between the barrier 
parameters and other junction characteristics — particularly RA, the 
product of junction resistance with area (a simple and extremely useful 


parameter) — are discussed. 


THEORETICAL tunnel conductance calculations have 
been carried out by Simmons [1]. These have been 
extended and discussed by others, in order to deal with 
the calculation of barrier parameters using the fitting 
coefficients derived from electron tunneling conductance 
curves. The parameters are the barrier thickness s, the 
average barrier height ¢ (¢ = 3[¢, + ¢2] at zero bias), 
and Ad = $2 — ¢, the difference of the barrier heights 
at the two electrodes. Using the methods of Brinkman, 
Dynes, and Rowell [2] (BDR), and Albrecht, Keller, 
and Thieme [3] (AKT) parameters have been calculated 
for a sample of approximately fifty junctions with 
various barriers (e.g. oxides, organics, and oxides plus 
organics). A comparison of these parameters, along with 
fitting methods and data handling are shown. Since 
independent measurements of s [4], @ [5] and Ag have 
not been made, little can be said about a comparison 
between experimental and calculated values. There do, 
however, exist relationships between the calculated s, @ 
and A@ and simple junction parameters, such as RA 
(where R is the junction resistance and A its area) which 
provides a quick characterization of the junction pro- 
perties. These will be discussed. 

Junctions used in this study were prepared by first 
evaporating an aluminum base film onto a clean glass 
substrate. Exposure of the base film to a plasma dis- 
charge of oxygen, methane, ethylene, fluoroform, or 
oxygen followed by methane formed the barriers. 
Evaporation of lead top electrodes completed the 
junctions (areas ranged from 0.06 to 0.72 mm? and 
resistances from 10 to 2000 92). Measurements were 
taken at 4.2 K using an updated version of the bridge- 
minicomputer system described by Adler and Straus [6]. 


* Supported in part by the National Research Council of 
Canada. 


To calculate barrier parameters, we used the coef- 
ficients by, b; and b» of a polynomial fit to the normal- 
ized tunneling conductance curve o(V)/o, =b) + b\V + 
b,V? +... where o(V) is the derivative of current with 
respect to voltage and oy = o(— 30 meV). This normal- 
izing bias was chosen because all junctions in this study 
had superconducting Pb cover electrodes. In order to 
avoid large structure due to Pb phonons and the energy 
gap, our sweep avoided the + 30 meV region. Because 
of this, R is taken as 1/oy,. The fitting region used was 
typically 150—500 meV (using a 6th degree polynomial 
for convenience). This avoided the zero bias region and 
any low energy inelastic structure, thus concentrating 
on the elastic part of the conductance [Fig. 1(a)]. The 
fitting of oxide or oxide plus methane data was fairly 
easy due to a lack of prominent inelastic tunneling 
structure above 150 meV [Fig. 1(a)]. Data from junc- 
tions with organic barriers such as ethylene presented a 
problem because of the large conductance step at 
360 meV (due to C—H stretching). This was solved by 
lowering data above the step by the step height and 
then performing the fit; two regions were used in this 
case 185—345 meV and 380-499 meV [Fig. 1(b)]. The 
same procedure was used for fluoroform and methane 
data. By using the minicomputer system, data handling 
and polynomial fits could be carried out quickly and 
easily. 

Once the coefficients were extracted, barrier para- 
meters were calculated using the methods of BDR and 
AKT. The difference in the numbers obtained by the 
two methods for s and ¢ is negligible. For example, the 
values of Sax which turn out to be larger than those 
for Sgpr deviate approximately 0.03 A over a range of 
10.5—14.5 A. For ¢ the maximum deviation is about 
0.002 V from ~ 2.0—3.5 V and ¢axr is again greater 
than ¢ppr. Since the earlier work of BDR gives 
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Fig. 1. (a) Polynomial fit for an Al—oxide—Pb junction. 
The fitting region is from 150—500 meV. (b) Polynomial 
fit for an Al—ethylene—Pb junction. The fitting regions 
are 185—345 and 380—5S00 meV to avoid the 360 meV 
conductance step. Beyond 380 meV the polynomial fit 
parallels the conductance curve. 


satisfactory parameters, the rest of the numbers used in 
this paper will be those calculated using this method 
[(BDR method I — WKB approximation using equation 
(7)]. 

When calculating a set of parameters for comparison 
or graphical purposes, one should use the same general 
fitting region for the whole data set. Depending on the 
fitting region, the parameter values change. For a 
particular fitting region, however, the same trends 
among the parameters remain. 

In Fig. 2, a plot of ¢ vs s for Al—ethylene—Pb data, 
fit in the 185—345 and 380—S00 meV regions is shown, 
Chosen for a wide range of s and ¢ values, the data 
points, at first glance seem scattered with a general trend 
of decreasing ¢ as s increases. Looking at RA values, 
however, one sees that RA increases from left to right 
along the arrow direction, with data points of the same 
RA value (e.g. 11.0—11.5 2-mm?) lying along a line 
perpendicular to the arrow direction. These lines are 
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2.0 


11.0 12.0 13.0 14.0 15.0 
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Fig. 2. This figure shows @¢, the average barrier height 
plotted against s the barrier thickness. The numbers 
beside the data points give their RA (Q nm’) values. 
Points with the same RA values define straight lines 
parallel to each other. Going from left to right along the 
direction normal to the lines (indicated by the arrow) 
one sees that RA increases. There is thus some order 

to a seemingly scattered plot. 


actually gentle curves described by so'’? = C with the 
constant C dependent on RA. How to determine the 
constant will be shown later. 

The exponentizl dependence of tunneling current 
on barrier thickness can be demonstrated by plotting 
s'’? against In (RA). As shown in Fig. 3 this yields a 
straight line which is fairly independent of the barrier’s 
chemical composition; data includes oxide, oxide plus 
methane, fluoroform, methane, ethylene, and also 
junctions with contamination (e.g. formic acid recog- 
nized by its spectrum) or instabilities [7]. When present, 
an instability shows up as a noisy patch in the da/dV 
curve. The instability region for the 0-500 meV range 
is below RA = 18. Also, by looking at this graph one 
can determine the RA dependent constant C, previously 
mentioned, such that s¢!/* = C. The equation for the 
line determined by the data (with the exception of the 
ethylene data which tend to deviate from it) is given by 
so'/? = 0.961 In (RA) + 17.33. For junctions with 
different base or cover electrodes data points do not 
generally fall on this line. 

The data also show a relation between Ad/¢ and 
Ao/&, where 

_ o(+ 500) o(— 500) 

Me On On é 
and 
_ 1]o(+ 500) ‘i o(— 500) 


Z On On 


using a normalized conductance curve. This can be 
thought of as a steepness difference in the conductance 
curve for two bias polarities and may also be written as 
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Fig. 3. A graph of s¢'/? vs In (RA). Data is shown for junctions with barriers formed in a plasma discharge of oxygen, 
methane (CH,), oxygen followed by methane, fluoroform (CHF;), or ethylene (C,H,). 
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Fig. 4. This figure shows the linear relationship between 
Ad/¢ and Ao/G. It should be noted that the ethylene 
points which show more scatter than the rest of the data 
are also those which deviate from the line in Fig. 3. 


Ao _ 2[o(+ 500) — o(— 500)] 
6  o(+ 500) + o(— 500) © 


The resultant linear plot is shown in Fig. 4. Looking at 
this graph, the ethylene points show noticeably more 
scatter than the rest of the data. This fact, along with 
the deviation on the s¢'’* vs In (RA) graph (the points 
which deviate in Fig. 3 also do so in Fig. 4) suggests that 
there is something different about ethylene. 


If one returns to Fig. 3 one sees that this is not the 
plot one would use specifically for plotting the data 
[ie. according to BDR one would use In (RA) vs 
1.025so'/? + In (s/@'’*)] but rather a simplified 
expression which works quite well for the majority of 
junctions. Plotting the complete expression, one 
obtains the expected straight line with no data point 
deviation. The term thus responsible for the ethylene 
deviation involves s/¢'’?. For the other junctions this 
factor is fairly constant, while for ethylene it tends to 
be larger. Comparing our expansion with the BDR 
expression for parameter calculation one finds that 
b> = B(s*/¢) (i.e. the quadratic coefficient is pro- 
portional to s?/¢) where B is a constant. For the 
deviant ethylene points b, is larger than for the other 
barrier types. Going back to Fig. 4 one sees that for the 
more scattered ethylene values, A@/@ is smaller than 
expected for a given Aa/6 value. Using b;, one can 
express Ad/¢ as Db, (¢'/?/s) where D is a constant. In 
this expression, the inverse of s/¢'’? is involved and thus 
Ad/¢ will be smaller than expected. If one uses a value 
of s/o’? which for other junctions is fairly constant to 
calculate Ad¢/¢, one finds that the scatter of Fig. 4 is 
reduced. 

The way that the larger quadratic coefficient, b>, 
manifests itself is in the shape of the conductance curve 
(in the fitting region). The shapes of ethylene curves are 
different from those of the other barrier types. Whereas 
the curves of the others exhibit a fairly smooth curved 
shape of conductance in the fitting region used, the 
ethylene curves do not show this type of behavior over 


> . 7 
pee ae he ae ei na ot 


oo 
* 
‘ & 
& 
‘ 
Pe 
“< 
j i f ioe er . : 
ahi 
4 t 
SOY Pa w71 fog’ FBhr An - arn. € yt cate £8 pie eT i 2 Gaks 


“te a alga laa ® qu esd i laaw She i 
+ eed Sat DEM ae Ey es Asphaad He 
baakliaeee thther ye te Bap ee My 
{ j Ty, dy DEAN a Gare OS 
) youl eters esa 4 gis eet, 
. Uy Ohl Leo a iin) "i 
+9" S12 at Big Mi : Lot ae & call aire a ME oy 
ctl Hibs ce hs yb Oe ny waite . } AEs. 
43 Ob As tials nw Hebe’ rete ate alee a? s co 
MPa Ahow ty Re aes (ae wn to? rupee 4 a a ‘ 
PT Di ey ah Lue jai Ps f “4 1 Sal JRE" Sons 
AD el fli. eh a vlgrsen +E Be ahah eth ne 
ig # oak tined hg a ied ih - wa 2 


Sa aa JO) ut (aye aS ie ye veils lonie 7 
i. I 


Fin ph Hee iat ah gihe Naf ay atid. aeeey.% 
aff a tate weak bate is res ‘eae sae! 
* gah ps ti ‘e aa lle Walks ane ty phpaien: "le 
dr Sin elem ae ean ‘Tis 
super ete FF seit. pee is ro y 
als uG ae i Wa sean: sna: sable " ad Iti b 
ov PAR 9 wae sty i fo hie oily 
an wl oe at atts Biaee alias hare 


i 
eT) Gis 


“a score is ain st aia ay 


—— corr calle ra a 
iF . apie id wp hy rm oe me 


3a r ‘ sails aed 
ery ae toners r — f a 


% 


‘ i 
; yy’ 


ar hurte't 
U é 


al a 


952 


the whole range. but have instead a flatter region be- 
tween 185 and 300 meV. Indeed, if data are fit in the 
380—S00 meV region where the smooth curvature does 
exist, the deviation in Fig. 3 is eliminated and the scatter 
is Fig. 4 reduced. The cause of the different curvature is 
not known, however it may be due to a difference in 
barriers formed from the plasma discharge products of a 
double-bonded organic (e.g. ethylene) as opposed to 
those produced from a single bonded organic material 
(e.g methane or fluoroform). It would be a project by 
itself to determine what is actually happening. 
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In conclusion, one finds that for a wide range of 
junction resistances and area, the conductance can be 
represented by: oy =kA e 8% where k =e'®™ and 
6 = 1.041. The constants k and £ are not related to the 
material forming the barrier but rather to the electrode 
materials. The scatter observed for ethylene barriers is 
due to a relatively large quadratic coefficient resulting 
in a value of s/o'/* which is larger than the fairly 
constant value obtained for other barrier types. Finally, 
a linear relation exists between Ad/¢ and Aa/a. 
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Annealing effects in tunnel junctions (thermal annealing) 
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This paper describes the results of subjecting Al-oxide-Pb tunnel junctions to thermal annealing. 
By looking at the junction resistance, barrier parameters (thickness, average barrier height, and 
the separation between the barrier heights at the Al and Pb electrodes), and inelastic electron 
tunneling (IET) peak intensities, one can study time-dependent effects occurring in the junction. 
These observations indicate changes in the tunneling barrier and IET spectra, and suggest that 
surface hydroxyl (-OH) groups experience reorientation. 


PACS numbers: 73.40.Gk, 74.50. + r, 73.40.Rw, 68.20. + t 


In this paper we set out to explore the nature of the 
insulating barrier in metal-insulator-metal tunnel junctions. 
It is known that tunnel junctions left at room temperature 
for any length of time usually undergo changes in resis- 
tance.’ In order to prevent this, the junctions are stored ata 
low temperature (i.e., liquid nitrogen at 77 K). To carry out 
useful studies of phenomena such as chemisorption of organ- 
ic molecules on oxide surfaces, it is necessary to have some 
knowledge of the temporal behavior of the oxide barrier in 
junctions prior to chemisorption. This is particularly appli- 
cable to cases in which one is studying chemical changes 
occurring within the barrier (e.g., chemical reactions of mol- 
ecules infused into oxide junctions’). 

This work describes what happens to Al-oxide-Pb tun- 
nel junctions during thermal annealing. By looking at junc- 
tion resistance, barrier parameters (thickness, average bar- 
rier height, and the difference between the barrier heights at 
the Al and Pb electrodes), and inelastic electron tunneling 
(IET) peak intensities, one can study time-dependent effects 
occurring in the junction. One observes changes in the tun- 
neling barrier and IET spectra indicating that surface hy- 
droxyl groups ( — OH) experience reorientation. 

The Al-oxide-Pb junctions used in this study were pre- 
pared in a vacuum system which was first cleaned with a ~ 5 
min oxygen discharge followed by ~ 15 min of argon plasma 
cleaning. An aluminum base film (2000 A thick) was then 
evaporated onto a clean glass substrate. Oxidation of the 
base film was carried out in an oxygen-helium (40-60%) 


FIG. 1. This figure shows a plot of R,/R, versus t for several junctions, 
where R, is the junction resistance at time #, Ro is the initial resistance, andt 
is the total annealing time. — 
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plasma [using Matheson research-grade oxygen (99.99% 
purity) and Helium (99.9999% purity)] with a total gas- 
mixture pressure of about 200 4 Hg. To create the plasma, 
two circular concentric copper electrodes were used, the out- 
er one being a circular band of copper sheet and the inner one 
a circular piece of copper rod. Typical plasma currents 
ranged from 2 to 6 mA with voltages of the order of 375-425 
V (outer electrode positive, inner electrode negative) for a 
discharge time of about 30 min. Two lead electrodes (2000 A 
thick) were then evaporated (at a rate of ~900 A /min.) to 
form two junctions on the same substrate sharing a common 
base layer (one of which could be used as a control junction 
when required). The entire fabrication was done without 
breaking the vacuum. Once formed, junctions were cooled to 
77 K in liquid nitrogen until ready for study. Measurements 
were taken at 4.2 K using an updated version of the bridge- 
minicomputer system described by Adler and Straus.’ 
Having briefly described junction fabrication, one 
should note that the insulating barrier formed in a plasma is 
not likely to be stoichiometric alumina and may contain oth- 
er groups (e.g., -OH). The covering metal electrode will per- 
turb the barrier in at least two ways: It may penetrate into the 
barrier, or it may alter the orientation of surface groups (i.e., 
sit on or squash them somewhat). Also, due to the work 
functions of the two-metal electrodes, a strong electric field 
will arise across the barrier* which will have some effect. For 
an idea of the field strength, a potential difference of 0.5 V 
(similar to differences in work functions) across a 10-A junc- 
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FIG. 2. A typical IET spectrum for an Al-oxide-Pb tunnel junction is 
shown. Peaks of interest occur at 118 meV and 450 meV. 
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_ FIG. 3. The change in the 450-meV peak with annealing is shown here for 
two different times: 0 and 32 h. 


tion gives a field of 5x 10° V/m. This may aid in group orien- 
tation and ionic diffusion. 

The thermal annealing experiments were carried out by 
removing a junction from liquid-nitrogen temperature and 
placing it at room temperature, about 20 °C, for the desired 
length of time. To prevent the ice which formed on the tun- 
nel junction electrodes (when removed from liquid nitrogen 
to room air) from melting and subsequently damaging the 
junction, annealing was carried out in a vacuum. All anneal- 
ing times referred to in this paper are cumulative times (i.e., 
total time the junction has been annealed). The most obvious 
parameter which changes with annealing is the junction re- 
sistance R. Figure 1 shows a plot of the resistance ratio R,/ 
Ro versus ¢ (total annealing time) for a series of thermally 
annealed (open circuit) oxide junctions where R, is the ini- 
tial resistance (within a few minutes of junction formation), 
and R, the resistance at time ¢ (at — 30 mV). The curves 
Suggest three regions: (1) An initial decrease (some curves 
drop below R,/R, = 1.00) in resistance (or none at all), (2) 
A rapid resistance increase which becomes, (3) a more gentle 
increase. 

Other noticeable changes occur in the conductance 7 
versus applied bias V, and in the second derivative 
do/dV = d*I /dV’, the IET spectrum. The spectrum of an 
oxide (Fig. 2) shows two main peaks of interest, the 450-meV 
peak due to surface OH stretching, and an Al-O stretching 
mode at 118 meV due to the aluminum oxide (all peak posi- 
tions referred to in this paper are for the polarities Al nega- 


FIG. 4. The model for the trapezoidal barrier with thickness s and barrier 
heights ¢, and ¢, (at the Al and Pb electrodes, respectively) is shown here, 
with an applied bias eV. 
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FIG. 5(a) A plot of the separation between barrier heights at the two elec- 
trodes 4¢ versus total annealing time ¢ shows the decrease of 4é with time. 
(b) This plot of the barrier heights at the two electrodes versus time shows a 
decrease at the Pb side (4, ) but little change at the Al side (¢, ). (c) This 
graph shows how the 450-meV peak intensity increases with annealing time 
and finally levels off. 


tive and Pb positive). Upon thermal annealing one often sees 
noticeable changes in the spectrum of a particular junction: 
For example (Fig. 3), the OH mode at 450 meV increases in 
size and sharpens up slightly. 

Some changes also occur in the shape of the conduc- 
tance curves. One can conveniently characterize these 
curves’ by calculating the barrier parameters, with s, the 
average barrier thickness, ¢, the average barrier height 
[¢ =4(¢, + ¢,) at zero bias], and Ad = ¢, — ¢,, the differ- 


thermal 


anneal 


FIG. 6(a) A schematic representation of the tunnel barrier shows how the 


‘OH groups reorient under the influence of an electric field during thermal 


annealing. (b) Here, a close-up view shows how @ decreases with reorienta- 
tion and results in the hydrogen side of the OH group moving closer to the 
Pb electrode. 
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FIG. 7(a) A plot of the 118-meV (AI-O) intensity versus annealing time for 
several tunnel junctions. (b) This graph of the 450-meV (-OH) intensity 
versus annealing time for several tunnel junctions shows that for all the 
samples the same saturation intensity is reached. 


ence of the barrier heights at the two electrodes (Fig. 4). This 
was done assuming a trapezoidal barrier and the WKB ap- 
proximation following the method of Brinkman, Dynes, and 
Rowell‘ (Eq. 7). The model is not taken literally but rather as 
an indicator of the qualitative changes taking place in the 
barrier. For thermally annealed junction, one finds that s, 
the effective barrier thickness, increases with time (this does 
not mean the barrier is necessarily becoming physically 
wider but rather the tunneling resistivity is becoming great- 
er), and that ¢, the average barrier height, decreases with 
time. The most interesting results involve the barrier heights 
at the Al and Pb electrodes, ¢, and ¢,. Figure 5(a) shows a 
decrease in Ad which levels off within 40 h. In Fig. 5(b) one 
sees that ¢,, the barrier height at the aluminum electrode, 
remains fairly constant whereas ¢,, for the lead electrode, 
decreases and then levels off. This suggests that there has 
been an increase in positive charge density near the lead elec- 
trode. Another interesting quantity to study is the area un- 
der the peaks of interest in Fig. 2 referred as peak intensities. 
It should be noted that our data is calibrated and we thus 
work with absolute values of da/dV = d7I/dV’, as op- 
posed to raw second-harmonic signals. From the data, the 
integrated peak intensity is calculated by 


Fv) = | Je, ar - «(| av, 


where a, is the conductance used to normalize do/dV [in 
this paper 7, = o(Vy) = o{ — 30 meV), since we did not 
sweep through zero bias in order to avoid structure due to 
the superconducting energy gap and phonons of Pb], and 
g(V) is the smooth background under the inelastic peak.’ 
One finds that the intensity of the 450-meV peak for the 
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PLEAS] 


junction described increases with time and tends to saturate 


[Fig. 5(c)]. 
Looking back over the results we find 


(1) A rapid 450-meV peak intensity increase to 
saturation, 

(2) A decrease in Ad, 

(3) A decrease in ¢, (the lead side becoming more 
positive), 

(4) d, remains constant for the aluminum electrode, 

(5) A slight sharpening of the 450-meV peak. 


Since there is a large electric field in the barrier ( ~ 10° 
V/m), a possible interpretation can be given by considering 
the reorientation of the surface hydroxyl groups on the alu- 
minum oxide. The above mechanism is shown schematically 
in Fig. 6(a) where the originally quite random OH groups 
reorient themselves under the influence of the electric field. 
Figure 6(b) shows that, in realigning the OH dipole, O (the 
average angle of an OH group from the normal to the plane 
of the Pb electrode) tends to 0°, i.e., the OH stretch motion is 
more perpendicular to the Pb electrode. This will result in an 
intensity increase for this mode proportional to cos” O for a 
stretching mode.*® Secondly, Fig. 6(b) shows that with reor- 
ientation the positive hydrogen of the OH moves closer to 
the lead electrode, increasing the positive charge density 
near it; therefore ¢, decreases while practically no effect is 
exerted on ¢, at the Al side. 

One also finds that upon annealing, the intensity of the 
118-meV peak due to Al-O stretching increases in a manner 
similar to the 450-meV increase. By looking at the intensity- 
vs.-annealing time curves [Fig. 7(a)] for a few junctions, one 
can see that the lower the initial 118-meV intensity the great- 
er the annealing effect (as measured by the ratio of saturation 
intensity at ~ 60 h to initial intensity), suggesting that more 
reorientation can take place the lower the initial 118-meV 
intensity. Whereas the 118-meV intensities tend to saturate 
at different values, the 450-meV intensities all tend to satu- 
rate at approximately the same value [Fig. (7(b)]. This sug- 
gests that there is a certain “monolayer” coverage of surface 
OH for all junctions (which has its own particular initial 
orientation) while the amount of aluminum oxide [Fig. 7(a)] 
can vary, 

To summarize, we have found that junctions subjected 
to thermal annealing undergo resistance changes (up and 


down) as well as changes in barrier parameters and IET peak 


intensities with time. This suggests reorientation of the insu- 
lating barrier (surface OH and aluminum oxide) under the 
influence of a natural electric field. Also, by looking at the 
values at which the IET peak intensities saturate, one de- 
duces that there is a certain monolayer coverage of surface 
OH for all the junctions while the amount of aluminum ox- 
ide varies. 
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Annealing effects in tunnel junctions (voltage annealing) 


M. K. Konkin and J. G. Adler 
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This paper describes the time-dependent changes which occur in Al-oxide—Pb tunnel junctions as 
a result of voltage annealing (i.e., keeping a fixed bias across junction electrodes). By examining 
changes in junction resistance, barrier parameters (barrier thickness, average barrier height, and 
the separation between barrier heights at the two electrodes) and inelastic electron tunneling 
(IET) peak intensities, one finds that applying a voltage across the junction causes an 
accumulation or depletion of positive charge at the Al-oxide interface depending on the polarity 
of the applied bias. Along with charge redistribution, the voltage-induced electric field also 
enhances or retards other processes (such as reorientation of surface hydroxyl groups) which 


normally take place during thermal annealing. 


PACS numbers: 73.40.Gk, 74.50. + r, 73.40.Rw, 68.20. + t 


Previous work’ has described some of the effects which 
occur in Al-oxide—Pb tunnel junctions when they are ther- 
mally annealed. The current paper describes the results of 
voltage annealing. In this, a constant voltage is applied 
across the junction electrodes to produce an electric field 
which adds or subtracts from the junction’s intrinsic electric 
field. This intrinsic field is responsible for some of the 
changes observed in thermal annealing; thus, altering the 
field with an applied bias should help to clarify our knowl- 
edge of the annealing process. Information about junction 
resistance, barrier parameters, and inelastic electron tunnel- 
ing (IET) peak intensities as they change under the influence 
of voltage annealing indicates that positive charge at the Al- 
oxide interface can be increased or decreased, depending on 
the polarity of the applied bias. Also, by proper choice of bias 
polarity, junction aging (i.e., any changes occurring within 
the junction after its completion) can be retarded. 

The Al-oxide—Pb tunnel junctions used in this study 
(about 20) were prepared by first evaporating an aluminum 
base layer onto a glass substrate. Oxidation of the aluminum 
film in an oxygen-helium plasma discharge’ to form the in- 
sulating barrier was then followed by the simultaneous evap- 
oration of two lead top electrodes. In this manner, one junc- 
tion of the pair served as a control (thermally annealed only) 
while the other junction was subjected to voltage annealing. 
Newly completed junctions had R,A values (proportional to 
junction thickness) ranging from 15 to 55 2 mm,’ where R, 
was the resistance of the virgin junction and A was the area. 
They were stored in liquid nitrogen, while measurements . 

’ were taken at 4.2 K using an advanced version of the bridge- 
minicomputer system described by Adler and Straus.” 

Having briefly described junction fabrication, one 
should note that the insulating barrier formed in the plasma 
is not likely to be stoichiometric alumina, and may contain 
other groups (e.g., — OH). The covering metal electrode 
will perturb the barrier in at least two ways: the electrode 

_ may penetrate into the barrier,’ or it may alter the orienta- 
tion of surface groups (e.g., surface OH) canting them some- 
what. Also, due to the work functions of the two metal elec- 
trodes, a strong electric field will arise across the barrier.* 

For an idea of field strength, a potential difference of 0.5 V 
(similar to differences in work functions) across a 10-A junc- 
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tion (these are typical barrier thicknesses found using ellip- 
sometry’) gives a field of 5 10° V/m. This field can aid in 
group reorientation and ionic diffusion. 

The voltage annealing experiments were carried out at 
room temperature in a vacuum. The control junction was 
allowed to undergo thermal annealing (open circuit) while 
the other junction of the pair was biased at the desired volt- 
age and polarity. We shall refer to a bias with A] positive, Pb 
negative as forward polarity (FP) and Al negative, Pb posi- 
tion as reverse polarity (RP). Previous thermal annealing _ 
experiments showed that the intrinsic field of a junction was 
of forward polarity. 

Interpretation of the annealing process is carried out by 
examining changes in junction resistance, barrier param- 
eters, and inelastic electron tunneling peak intensities with 
respect to total annealing time. The junction resistance R, is 
the measured resistance (at —30 mV) at time ¢ (total anneal- 
ing time). Changes occurring in the shape of the conduc- 
tance curves (o versus bias V for a sweep range of + 500 
mV) can be characterized® by calculating the barrier param- 
eters s, the average barrier thickness ¢, the average height 
[ d = 1(¢, + ¢) at zero bias], and Ad = ¢, — ¢,, the differ- 
ence of the barrier heights at the two electrodes (values for 
the parameters ¢, and ¢, can be easily determined from the 
already calculated ¢ and A¢ ). This was done assuming a 
trapezoidal barrier and the WKB approximation following 
the method of Brinkman, Dynes, and Rowell’ [Eq. (7)]. The 
model is not taken literally but rather as an indicator of the 
qualitative changes taking place within the barrier. Finally, 
the last quantities of interest in this study are the areas under 
peaks of interest in the IET spectrum do/dV referred to as 
peak intensities. From our calibrated data, the integrated 
peak intensity J is calculated by 


r= { |2-% —av) |av, 
Gy. dV 

where cy is the conductance used to normalize da/dV [in 
this paper a, = o(V,,) = of —30 meV) = 1/R,, since we 
did not sweep through zero bias in order to avoid structure 
due to the superconducting energy gap and phonon strurc- 
ture of Pb] and g(V ) is the smooth background under the 
inelastic peak.? 
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FIG. 1. This figure shows a plot of R,/R, vst — t, for three junctions biased 
at 0.5-V FP. R, is the junction resistance at time t, R, is the resistance at the 
time of bias application ¢,, and ¢ is the total annealing time. For the circles, 
squares, and triangles t, = 0, 212, and 324 h, respectively. Each graph point 
is the value measured after several hours of 0.5-V FP annealing. 


A number of observations were made in the course of 
the annealing studies. Changes in the control junction resis- 
tance, barrier parametes, and peak intensities were consis- 
tent with previous results.’ In particular, reorientation of 
surface OH groups as monitored by the 450-meV peak (due 
to surface OH stretching) was suggested by (1) peak intensi- 
ty increase to saturation, (2) a decrease in Ad, (3) a decrease 
in ¢, (the lead side becoming more positive), and (4) little 
change in ¢, for the aluminum electrode, all with respect to 
time. This reorientation of surface OH as well as the alumi- 
num oxide (indicated by Al-O stretching at ~ 118 meV) oc- 
curred due to the intrinsic field of the junction trying to align 
the dipoles. Voltage annealed junctions showed several new 
effects. 

Junctions biased at 0.5 V with forward polarity (Al 
positive, Pb negative) show an increase in the slope of the 
resistance versus annealing time curves shortly after the de- 
sired voltage is applied. In Fig. 1 a graph of R,/R, versus 
time t — t, (where R, is the junction resistance at timet, R, is 
the resistance at the time of bias application t,, and ¢ is the 


0 25 50 75 


FIG. 2. A plot of 4g — 4¢, vst — t, for three junctions shows the effect ofa 
0.5-V FP bias on A¢, the difference of the barrier heights at the two elec- 
trodes. The junction age when biased t, equals 0, 212, and 324 h for circles, 


squares and triangles, respectively. A¢, is the value of A¢ at time ¢,. 
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FIG. 3. The resistance behavior R,/R, vs time is shown here for a junction 
biased at 0.5-V RP (circles), and a thermally annealed control junction 
(squares). 


total annealing time) illustrates this resistance increase for 
several junctions of different ages when initially biased. One 
sees that the increase soon slows down again and that the net 
resistance ratio increase after, say 50 h of application of the 
bias, diminishes the older the junction was to start with. 

Another noticeable change is that which occurs in ¢, 
and ¢,—the potentials at the Al and Pb electrodes, respec- 
tively (depending on the nature of the sample ¢, can be great- 
er or smaller than ¢,). Following application of an annealing 
bias, it is found that Ag = ¢, — ¢, increases in time and fi- 
nally levels off as shown in Fig. 2 by Ag — Ad, vst —t, 
(4¢, refers to the value of A4¢ at time ¢, ). It is seen also that 
an age dependence is exhibited by the magnitude of the ef- 
fect. One should note that it is 6, which decreases markedly 
while ¢, changes only a little. This implies that the barrier at 
the aluminum electrode acquires positive charge with time 
which eventually saturates as long as the biasing voltage is 
maintained. This charge, moreover, is trapped just within 
the barrier otherwise, like a capacitor, it would drain off as 
the annealing voltage was reduced to zero. 

If junctions are annealed using reverse polarity (Al neg- 
ative, Pb positive) at 0.5 V, the results are different. Looking 
at R,/Ry vst (Rois the resistance value at t = 0), one finds an 
initial decrease in resistance which eventually increases 
slowly (Fig. 3). The barrier height at the Al electrode ¢, goes 
up (implying that the Al side has become less positive) while 
that at the Pb electrode ¢, remains constant. 

These resistance changes as well as the behavior of ¢, 
(Fig. 4) and ¢, for a biased junction suggest movement of 
positive Al ions in and out of the barrier at the Al-oxide 
interface. When forward polarity is utilized, aluminum ions 
may be created at the aluminum electrode and drift a short 
way into the oxide barrier under the influence of the applied 
field. Such an occurrence raises the amount of positive 
charge near the Al electrode (¢, thus decreases) and in- 
creases junction resistance since the inclusion of Al ions 
within the oxide makes the barrier denser near the Al-oxide 
interface. On the other hand, application of a reverse polar- 
ity bias causes an initial resistance decrease and an increase 
in ¢,. This indicates that some already existing positive 


M. K. Konkin and J. G. Adler 5451 


Aon 
ie 
| ‘". 
‘si 
| a. 
¥ 
- “vv 
a, 
Fiend eer tad rivet sane 
sieseoaanieit acts 4 1 elie 
‘S16 Sacre elant, port yt peer it gerbe 
(oni). sateay| nae -pewinempery 19S bas LB 
qtttenein «hi wmind 9 
‘ :) fans outlets lei 7) 
New SST 1S ect Lae eat fay ai ath ah! Pina 
th COG Up elon oid Pap ee ier oat i atneiey 
Oe Vik FIAT (hier tapas | 
a i: nm! aiicg od? Tole Rae 
f h fis *) 
a , Eb rice al tacos 
i “ tie ae ev Art ea? 
ye? es eG lid tad 
» Mi ed fh aera By aes Sen 
bein? 2s aly BE aie meee) caer, 
is of Tek accra) et LSCONies oF AStiaten gi ih a 
Stainien. mee nioy dn oA Sport Deets tne Sart toed 
ti oti jets Brel ale oe viet way, Me 3 ty Hit 7. 
weit igs Pate Huon culate aboees wil ees 
aor yy aries iW ad fred. te ev eres i peiias afr. fi tak, , +e : 
Py TEL dea, ies? ci Le Oma As; re itt aye aa Sree rare" 
ad “io orxd be 4 Meo s a peoini i, HOY 1 1 
aT (MoT eit RIN alae ch et. 
? vy -1z ee enh Tinea oboe a * f 
snuloa T veairth ber eelen state. e 013 tei eal Ae 
enol sig dG | tones sorenigns “tag: /2EV RS ihe 
sar ker Ves ees ica MEY ESI Nyy up Apis 


authe’, “> Abarercale Uvedls uit hgled “preaiedl tt ay ih 
etitey Covi’ gab sovoce tl duet able wetrigty 
Saadaiieese heen * arias ae 
>}: wwiteded eth og lof aa bo gndele pe Pas 
ERAT Upgaar apron bent? ro 
same eh aee stave cH hers bik oi ene Ai 
ann musioaglé Eastin at mon. Hite. pb ta 
yenke ns fipb.kby storia: ts au lianas sd? hea gi oon 
besblg te Ate 1G sepantthat a ‘Sit tirectatine Pcl 
SH, Me yagaers cd whe lar tes se malbed es sig a 
~tid ius Cagenancaly aud 2h) ae eal 2, anid vedi 
Moll Tq ondnbatl tt) seatbeuriiinns ental! res 
obit A als mean Brett c as) bet Salt wasn ne 
valg W999 Ae are 8 aye 
we onl; on bade 
avesilnng > er a 


egg 


ras StS iy | 
jan 


my 
“8. _8___8 


-—--g-----8---#s 


“0 20 4 60 80 100 


0 
t(hr) 


FIG. 4. This figure illustrates the behavior of the barrier height at the Al 
electrode ¢, vs time for junctions subjected to 0.5-V RP (circles), 0.5-V FP 
(squares), or thermal annealing only (triangles). 


charge is moved out of the barrier, which is reasonable since 
a negative potential has now been applied at the Al side. 
One can account for the age dependence of the amount 
of change in R,/R, and Ad by noting the behavior of a volt- 
age annealed (FP) junction once the annealing bias is re- 
moved. Figure 5 shows what happens. A decrase in Ad [Fig. 
5(a)] occurs and ¢, increases suggesting that the amount of 
positive charge near the Al electrode is being depleted. 
Looking now at the junction’s resistance curve [Fig. 5(b)], it 
is seen that as charge depletion occurs, resistance increases. 
This leads one to believe that Al ions are being neutralized, 
and remain within the barrier as opposed to drifting back out 
to the Al electrode. The age effect can now be explained as 
follows. As a junction anneals, Al ions are created and drift 
into the junction due to its intrinsic field or under the influ- 
erice of an externally applied field. The ions, once having 


0 50 100 


0 50 100 
t(hr) 


FIG. 5(a). Here we see the decrease in A¢ after the removal of a 0.5-V FP 
bias. The arrow indicates the beginning of thermal annealing. (b) The resis- 
tance behavior of a junction during voltage annealing (0.5-V FP) and then 
thermal annealing (starting at the arrow) is shown. 


5452 J. Appl. Phys., Vol. 51, No. 10, October 1980 


. 


| oxide Pb 


FIG. 6(a). This figure shows the response of a “young” junction to FP 
biasing. In the first picture many vacancies exist at the Al-oxide interface. 
When a bias is applied (second picture) ions are created and can drift into 
the barrier to fill the vacancies. They remain trapped here after the bias is 
removed (picture three). (b) Here, the response of an “‘old’’ junction to FP 
biasing is shown. Vacancies are filled the same way as described in Fig. 6(a) 
except that there are fewer vacancies to begin with. Circles with + signs 
represent Al ions, and shaded-in circles denote Al atoms. 


entered the barrier are eventually neutralized and new ions 
can be formed (maintaining a type of dynamic charge equi- 
librium). As a junction ages, vacancies for new ions to occu- 
py within the insulator are diminished since some have al- 
ready been filled. Thus, when an external bias is applied to an 
old junction as opposed to a newer one, its response is not as 
great simply because there are not as many vacancies in the 
barrier for new AI ions to fill and produce their effects. This 
is illustrated in Fig. 6. On the subject of aging, one finds that 
junction aging can be suppressed by the application of a re- 
verse polarity bias. We will discuss this later. 

Now, we turn to see what effects voltage annealing has 
on IET peak intensities. Looking at an oxide spectrum (Fig. 
7), one notes peaks of interest at ~ 118 and 450 meV due to 


Al-oxide-Pb 


-118 Al-O stretch 
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E (meV) 
FIG. 7. The two curves are IET spectra for an Al-oxide—Pb junction at two 
different times. The upper curve shows the effect that 70 h of 0.5-V FP 


voltage annealing has on the virgin sample (lower curve). Peaks of interest 
occur at 118 and 450 meV. 
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FIG. 8. This plot of peak intensity J vs time shows how reverse polarity (0.5 
V) annealing suppresses any changes in J (0-120 h). If the bias is removed 
(at the arrow) the junction thermally anneals and peak intensities increases. 


Al-O and O-H stretching, respectively. When voltage an- 
nealing is carried out with a 0.5-V FP bias, one finds the 118- 
and 450-meV peak intensities increase with time to values 
greater than those in the thermal annealing case. This sug- 
gests that with an electric field greater than the intrinsic field 
alone, the surface hydroxy] and aluminum oxide may exper- 
ience further reorientation. On the other hand, the applica- 
tion of a field of 0.5-V RP is found to inhibit the change in 
peak intensities which normally occurs in virgin junctions 
under the influence of their intrinsic field (Fig. 8). This im- 
plies that the intrinsic field of the junction is neutralized by 
the 0.5-V RP bias, thereby leaving no electric field across the 
junction to reorient the molecular groups within. If the bias 
Is now removed, the junction thermally anneals as before; 
resistance and peak intensities increase. By applying such a 
reverse field to neutralize the intrinsic field, one can thus 
delay the aging of a tunnel junction. 

Having looked at the voltage annealing behavior, one 
can now use it to aid in explaining some results observed in 
the thermal annealing experiments. A number of thermally 
annealed junctions initially showed an unexplained drop in 
junction resistance similar to the one observed in voltage 
annealing with reverse polarity. Unlike RP anneaing the 
peak intensities increased in the usual way. A close look now 
at ¢, for those junctions with a pronounced resistance dip 
shows that ¢, initially increases a little (Al side becoming 
more negative). This, as well as the resistance drop, could be 
accounted for if positive Al ions left the barrier. However, 

- since the intrinsic field of the junction has the Al electrode 
positive, there seems to be a contradiction unless one consid- 
ers what may occur in the course of junction fabrication. 
Looking at the barrier during growth, one might imagine it 
to be aluminum oxide with positive Al ions near the alumi- 
num electrode and OH groups standing up on the growing 
alumina surface [Fig. 9(a)]. Subsequent evaporation of the 
Pb cover electrode will tend to squash the OH groups, leav- 
ing them canted as shown in Fig. 9(b). This movement of the 
OH group’s hydrogen closer to the oxide will cause a greater 
electrostatic repulsion between the aluminum ions and the 
hydroxyl hydrogens than existed prior to the deposition of 
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the Pb cover electrode. Since the Al ions are smaller and 
more mobile than the OH groups, the mutual electrostatic 
repulsion can force some ions out of the barrier [Fig. 9(c)]. 
This takes place only until the intrinsic field of the junction 
can reorient the OH groups sufficiently to relieve the repul- 
sive force. At this point, Al ions under the influence of the 
intrinsic field are once more able to enter the barrier [Fig. 
9(d)] and cause resistance increases. Some thermally an- 
nealed junctions show the above behavior; some do not. Fac- 
tors influencing this probably include (1) the manner in 
which the Pb electrode is evaporated and (2) the number and 
distribution of Al ions within the oxide. 

In this paper, we have introduced a new technique for 
the study of ion motion and behavior within tunnel junc- 
tions. By applying a 0.5-V FP annealing voltage to a junc- 
tion, one causes positive charges, presumably aluminum 
ions to accumulate at the Al-oxide interface. These ions 
serve to increase Ag (¢, decreases) as well as the junction 
resistance. When the voltage is removed, neutralization of 
ions within the barrier is reflected by the fall of 4é (¢, now 
increases). The accumulated charge also serves to enhance 
the junction’s intrinsic electric field which further aids group 
reorientation. When a 0.5-V RP bias is applied, some posi- 
tive charge is removed from the barrier and junction resis- 
tance initially falls. It rises slowly afterwards, perhaps due to 
neutralization of any remaining charge within the barrier. 
Peak intensities do not change because the magnitude of the 
bias seems sufficient to negate the junction’s intrinsic field, 
thereby leaving no electric field to reorient molecular 
groups. Junction aging is thus postponed until the reverse 
polarity bias is removed, at which point normal thermal an- 
nealing occurs. Finally, an explanation for the resistance dip 
observed in some thermally annealed junctions has been sug- 
gested as due to a temporary electrostatic repulsion occur- 
ring between Al ions and the H of surface hydroxyl] groups 
after the Pb cover electrode has been evaporated. 


Sie telat 


e Eiieldo a 


FIG. 9. This figure shows how evaporation of the Pb electrode on the grow- 
ing oxide can perturb the surface hydroxyl groups and cause expulsion of Al 
ions out of the barrier due to elerctrostatic repulsion. 
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